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Abstract  of  “Steady  Waves  in  a  Nonlinear  Theory  of  Viscoelasticity 
by  Gregory  Thomas  Warhola,  Ph.D.,  Brown  University,  May  1988. 


This  work  considers  the  propagation  of  steady  waves  in  viscoelastic  materials 
for  which  the  nonlinear  strain  measure  is  not  necessarily  convex.  The  shape  of  such 
a  wave  is  governed  by  an  ordinary  nonlinear  integro-differential  equation  having 
a  possibly  singular  difference  kernel.  The  existence  and  structure  of  a  solution 
depends  upon  the  relation  of  the  wavespeed,  a  parameter  in  the  problem,  to  two 
speeds  based  upon  the  state  of  the  material  ahead  of  the  wave.  Solutions  are 
constructed  by  a  monotone  iterative  scheme  which  is  proven  to  converge  to  a  unique 
solution  within  restricted  classes  of  functions  depending  upon  the  wavespeed.  A 
simple  numerical  approximation  to  the  iterative  scheme  is  used  to  produce  graphs 
of  solutions.  An  algebraic  “quasielastic”  approximation  produces  upper  bounds 
on  discontinuous  (shock  and  acceleration  wave)  solutions.  For  a  material  such  as 
polymethyl  methacrylate  (PMMA)  having  a  small  power  in  a  power-law  model  of 
its  compliance,  this  approximation  is  found  to  be  useful  for  accurately  predicting 
the  structure  of  shock  solutions. 
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CHAPTER  1:  INTRODUCTION 
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We  are  concerned  with  the  propagation  of  mechanical  disturbances  in  visco¬ 
elastic  materials.  Any  such  a  disturbance  in  a  material  is  called  a  wave  if  it  has 
features  which  are  clearly  recognizable  as  it  progresses  through  the  material.  In  a 
one-dimensional  theory  where  x  and  t  represent  place  and  time,  respectively,  a  wave 
is  steady  if  its  features  are  described  by  a  single  shape  function,  say  w(9),  where 
6  =  t  —  x/V  .  The  constant  V  is  the  wavespeed,  the  velocity  of  propagation  of  the 
disturbance  to. 

Transient  loading  experiments  conducted  by  Schuler  [1.1]  on  the  viscoelastic 
material  polymethyl  methacrylate  (PMMA)  have  produced  waves  which  appear 
steady  over  intervals  of  observation.  In  this  and  other  viscoelastic  materials  having 
a  nonlinear  dependence  of  the  stress  upon  the  strain,  the  tendency  for  disturbances 
to  dissipate  due  to  time-dependent  viscoelastic  effects  is  opposed  by  the  nonlin¬ 
earity.  Shock  waves  are  produced  if  the  loading  is  sufficiently  high.  Such  waves 
are  mathematically  characterized  by  a  discontinuity  in  the  material  particle  veloc¬ 
ity.  Kolsky  has  produced  travelling  waves  in  stretched  natural  rubber  [1.2],  which 
waves  exhibit  a  rapid  variation  in  the  particle  velocity  at  the  wave  front. 


Early  mathematical  investigations  into  the  propagation  of  steady  acceleration 


and  shock  waves  were  conducted  by  Coleman,  Gurtin,  and  Herrera  [1.3],  in  which 


they  assumed  e<istence  of  such  solutions.  Acceleration  waves  are  characterized  by 


a  discontinuity  in  some  derivative  of  the  particle  velocity.  Such  waves  travel  at  a 


critical  value  of  wavespeed,  with  respect  to  the  state  of  the  material  ahead  of  the 


wave;  shocks  travel  at  supercritical  speeds.  Pipkin  [1.4]  demonstrated  the  existence 


of  steady  shock  and  acceleration  wave  solutions  for  a  specific  model  of  a  nonlinear 
viscoelastic  fluid.  Furthermore,  his  treatment  includes  smooth  solutions  at  sub- 
critical  wavespeeds.  In  an  abstract  mathematical  setting,  Greenberg  proved  the 
existence  of  steady  shock  waves  for  a  broad  class  of  nonlinear  viscoelastic  materi¬ 
als  [1.5].  Greenberg  and  Hastings  ji.6]  later  studied  steady  waves  in  viscoelastic 
materials  in  a  somewhat  less  abstract  setting.  A  lengthy  description  of  experimen¬ 
tal  and  mathematical  investigations  into  this  subject  until  1974  is  contained  in  the 
work  by  Nunziato  et  al  [  1.7]. 

Our  goal  is  to  study  the  propagation  of  steady  waves  in  nonlinear  viscoelastic 
materials  with  yet  less  abstraction  while  providing  a  treatment  which  is  sufficiently 
general  to  include  useful  current  models  of  such  materials.  This  is  a  purely  me¬ 
chanical  treatment;  no  thermodynamic  quantities  are  considered.  In  Chapter  2,  we 
provide  foundation  material  for  the  rest  of  this  work.  The  equations  governing 
material  deformation  are  presented  along  with  the  constitutive  equation  describing 
the  viscoelastic  materials  considered.  We  describe  at  some  length  the  general  non¬ 
linearity  we  consider,  since  it  includes  a  departure  from  the  convexity  requirement 
imposed  in  the  earlier  treatments  [1.5]  and  [1.6].  We  also  introduce  the  quasi-elastic 
approximation  to  the  constitutive  equation.  We  later  find  this  approximation  to  be 
accurate  enough  for  description  of  shock  waves  in  materials  like  PMM  A. 

We  derive  a  nonlinear  integro-differential  equation  governing  steady  waves  in 
Chapter  3.  This  equation  contains  the  wavespeed  as  a  parameter.  After  a  brief 
look  at  steady  waves  in  a  linear  viscoelastic  theory,  we  prove  some  general  results 
concerning  steady  waves  in  a  nonlinear  viscoelastic  material.  From  these  results,  we 
propose  candidate  solutions  to  the  problem,  which  solutions  depend  upon  the  rela¬ 
tion  of  the  wavespeed  to  a  critical  value.  In  Chapter  4,  we  prove  the  existence  and 


uniqueness  of  these  steady  wave  solutions.  We  use  an  iteration  scheme  to  produce 
monotone  sequences  of  successive  approximations  to  the  solution.  Our  approach  is 
in  the  spirit  of  the  work  of  Greenberg  and  Hastings  [1.5];  however,  we  have  relaxed 
some  of  their  hypotheses  on  the  materials  considered.  Additionally,  we  use  a  dif¬ 
ferent  iteration  scheme  which  is  amenable  to  a  simple  numerical  approximation  for 
the  construction  of  graphs  of  approximate  solutions. 

In  Chapter  5,  we  consider  the  detailed  structure  of  steady  shock  and  acceler¬ 
ation  waves  near  the  discontinuity.  We  also  discuss  the  continuity  of  such  solutions 
as  a  function  of  the  wavespeed.  We  show  how  the  information  in  a  known  solution 
at  a  given  wavespeed  can  be  used  to  obtain  the  solution  for  another  wavespeed. 
The  numerical  approximation  to  the  iteration  scheme  is  introduced.  We  use  it  to 
construct  approximate  solutions  for  shock  and  acceleration  waves  in  materials  char¬ 
acterized  by  exponential  and  power-law  moduli.  It  is  here  where  we  see  that  the 
quasi-elastic  approximation  is  close  to  the  exact  solution  for  power-law  materials 
having  a  small  power.  The  graphs  of  shocks  we  produce  are  in  good  agreement  with 
the  experimental  results  for  PMMA  obtained  by  Schuler  [1.1]. 

The  final  chapter,  6,  is  devoted  to  smooth  solutions  below  the  critical  (accel¬ 
eration  wave)  wavespeed.  We  use  a  perturbation  technique  to  obtain  solutions  for 
waves  of  infinitesimal  amplitude  in  a  viscoelastic  solid.  We  construct  solutions  via 
iteration  for  specific  examples  involving  power-law  materials. 


-  - 


Ms 


CHAPTER  2:  NONLINEAR  VISCOELASTIC  DEFORMATIONS 


2.1  Governing  Equations. 

We  consider  one-dimensional  deformations  of  an  infinite  homogeneous  nonlin- 
early  viscoelastic  body  from  a  (possibly  strained)  reference  configuration  in  which 
the  body  is  at  rest  and  its  material  points  (or  particles)  are  identified  with  their 
position  x  on  the  real  line.  Let  u(x,t)  be  the  displacement  at  x  at  the  time  t  of 
a  particle  from  its  reference  position.  We  consider  functions  u  which  are  continu¬ 
ous  and  piecewise  differentiable  in  each  of  their  arguments.  We  write  the  partial 
derivatives  du(: c ,  •)  /  d  x  =  ux(z,)  and  du(  ■,  t)  /  dt  =  uf(-,l).  The  perturbed  particle 
velocity  is 

v{x,  t )  =  ut(x,  t).  (2.1.1) 

The  perturbed  strain  is  given  by 

e(x,t)  =  ux(x,t).  (2.1.2) 

We  normalize  the  extra  stress  a(x.t)  by  the  mass  density  and  write  the  balance  of 
linear  momentum  as 

vt[x,t)  =  <7z(x,t),  (2.1.3) 

where  these  derivatives  exist.  From  equations  (2.1.1)  and  (2.1.2),  the  compatibility 
condition  is 

et(x,  t)  =  vz(x,  t).  (2.1.4) 

The  following  constitutive  relation  defines  the  nonlinear  viscoelastic  materials  con¬ 
sidered  in  this  work;  it  is  a  viscoelastic  extension  of  the  relation  a  =  f  (e)  used  in 


nonlinear  elasticity.  We  consider  materials  for  which  the  stress  at  the  current  time  t 
depends  on  the  current  strain  and  the  entire  previous  strain  history  through  the 
convolution  relation 

/+  00 

G(t  -  t  )df(e(x,  t)).  (2.1.5) 

—  OO 

The  integral  in  equation  (2.1.5)  is  a  Stieltjes  integral.  We  describe  its  properties 
which  are  essential  for  this  work  in  the  following  section.  All  of  the  nonlinearity 
in  the  problem  is  embodied  in  the  strain  curve  /(e),  whose  graph  is  local  to  the 
reference  state  such  that  /( 0)  =  0.  Further  properties  of  the  curves  /  that  we 
consider  are  contained  in  section  2.3.  The  function  G(t)  is  the  stress  relaxation 
function  with  units  of  modulus  divided  by  mass  density.  Its  properties  are  contained 
in  section  2.5,  after  a  section  on  regularly-varying  functions.  We  discuss  some 
properties  of  the  constitutive  law  in  section  2.6. 

We  note  that  in  a  more  general  treatment  one  could  consider  the  constitutive 
model 

/+  00 

[(G(l-  r)-  Ge)]d/(E(r,r)),  (2.1.6) 

-  OO 

where  Ge  is  a  positive  constant.  In  this  model  the  equilibrium  elastic  stress  Geg(e) 
has  a  different  nonlinear  dependence  on  the  strain  than  does  the  transient  integral 
term.  In  cases  where  the  curves  /  and  g  are  determined  experimentally,  this  ap¬ 
proach  doubles  the  data  required  to  model  the  nonlinearity  in  the  problem.  For 
this  work,  we  assume  that  the  nonlinearity  is  characterized  by  a  single  set  of  these 
measurements,  such  that  /  =  g,  in  which  case  equation  (2.1.6)  reduces  to  (2.1.5). 


1 


2.2  Properties  of  the  Stieltjea  Integral. 


Throughout  this  work  we  will  be  concerned  only  with  functions  which  may 
be  written  as  the  indefinite  integral  of  their  derivative,  in  the  generalized  sense 
described  below.  Such  functions  have  bounded  variation  on  any  finite  interval.  These 
properties  are  enough  to  ensure  the  existence  of  the  Stieltjes  integral.  A  proof  of 
its  existence  and  a  collection  of  its  properties  is  contained  in  the  monograph  by 
W  idder  [2.1],  We  illustrate  only  what  is  needed  herein.  Consider,  for  example,  the 


im  proper  integral: 


K  K  )<*p(t), 


(2.2.1) 


where  /j  is  monotone.  (If  p  is  not  monotone,  it  can  always  be  decomposed  as 
/i  =  /i  +  —  fi  _  ,  where  /i+  and  n-  are  monotone.)  If  n(t)  has  a  piecewise  continuous 
derivative  /z  7(< ) ,  then  dfi(t)  =  n'(t)dt  in  equation  (2  2.1),  which  is  then  evaluated 
as  an  ordinary  Riemann  integral.  Thus,  if  ft(t)  =  t,  then  /  reduces  to  the  Riemann 
integral  of  K  .  If,  however,  n  has  a  jump  at  t0  given  by 


A fi(t0)  =  lim  /z (<)  -  lim  /z(f), 
t  1  to  t\t„ 


/+  00 


(2.2.2) 


(2.2.4) 


In  particular,  for  a  function  J  (t)  which  vanishes  for  t  <  0,  has  a  finite  value  J0 
at  t  —  0,  and  whose  derivative  is  given  for  t  >  0  by  J  '+(t),  we  have 


dJ(T) 


=  /  dJ(T) 


—  j  0  + 


/  J+(r)dr 


=  J0  +  J  (t)  —  Jc 


=  din- 


(2.2.5) 


Integrals  like  /(|  dJ  (r)  are  understood  to  include  integration  of  the  jump  in  J  at  an 
endpoint  of  the  in  teg  rat  ion  interval,  if  such  a  jump  exists  We  use  the  not  a  t  ion  /  „  ( / 1 
to  denote  the  restriction  of  /  to  positive  values  of  t 

2.3  The  Strain  Curve  /(e). 

The  strain  curve  /  (e)  is  assumed  to  be  a  single -valued  continuous  function  of  - 
with  a  pier?  w  ise  con  tin  nous  derivative  f'(e)  When  the  deform  ations  are  sin  all  w  .• 
expect  that  the  linear  theory  of  viscoelasticity  would  apply;  this  leads  us  to  consider 
curves  /  wh  ich  have  the  property 

/(e)  -  e,  e  -  0  (2  3  l  ) 

We  limit  discussion  to  functions  /  which  are  strictly  greater  than  e,  over  the  range 
of  e  considered,  when  the  strain  is  non-zero: 

/(e)  >  e,  e  *  0 

/(e)  =0,  e  =  0  (23  2) 

We  do  not  require  /  to  be  globally  convex  Instead,  we  consider  the  possibly  non 
convex  functions  whose  graphs  have  a  unique  intersection  with  any  ray  from  the 
origin,  over  the  range  of  e  considered  For  the  cases  in  which  we  are  most  inter¬ 
ested,  €  is  positive.  We  require  that  the  equation 

f  (e)  -  m  e  =  0 ,  m  >  1  .  (2  3  31 

have  a  unique  positive  solution  c.  for  each  value  of  m  W' e  now  show  that  these 
properties  are  sufficient  for  us  to  write  any  /  we  consider,  for  positive  strain,  as 


/(c)  =  m  (c)  c  , 


c  >  0. 


Vi  3  4) 


where  m  |  f )  is  a  continuous  function,  increasing  m  onotoiiu  all>  from  m(u)  -  I 
vs  ■  l  li  a  |>  ice  ew  ise  continuous  derivative  in’ft)  We  already  li  .1  v  .■  as  an  assump- 
tion  (2  3  ;t  |  that  -  is  uniquely  determined  by  m  Non  tn  !\f)h  uniquely 

determines  m  as  a  function  of  f  for  t  >  0  Furthermore.  m|n|  1  follows  directly 

from  (2  3  I  |  Thus  there  is  a  one-to-one  re  lat  ion  sh  i  |i  between  in  and  <  m  ( e  I  is 
monotone  Property  (2  3  2)  implies  that  m  is  increasing  |strn  tly  since  rn  is  one-to- 
one)  T  h  e  con  t  ill  u  it  >  of  /  requires  that  m|t)  tie  rootinuiius  Since  »e  require  /  '|f  I 
to  I.,-  pi.  cewisc  continuous,  and  since 

/  '(c  )  =  m  (t )  +  m  '(f  )  f  (2  3  5) 

m,(-|  must  too  be  pierew  ise  continuous 

Conversely,  any  function  /  of  the  form  (2.3.4)  necessarily  has  the  properties 
we  desire  The  continuity  of  /  and  the  piecewise  continuity  of  its  derivative  follow 
from  the  like  properties  of  rn  and  m'  when  used  in  equations  (2  3  4)  and  (2  3  5) 
Properties  (2  3  1)  and  (2  3  2)  are  satisfied  since  m(0)  =  1  and  m  is  increasing 
Finally,  (2  3  3)  has  a  unique  solution  e  since  //e  =  m(e)  is  monotone 

A  representative  curve  /(f)  is  shown  in  Figure  2  1  From  (2  3  4).  we  have  that 
the  slope  of  the  secant  from  the  origin  to  the  point  (c,  /(e))  is  m(f)  Since  m(e)  is 
increasing  we  see  from  (2  3  5)  that 

/  '(e )  >  m  (e).  f  >  0,  (2  3  6) 

with  equality  if  and  only  if  m'(f)  =  0;  the  points  having  this  property  cannot  be  a 
dense  set  sine**  m  is  strictly  increasing 

Fquation  (2  3  6)  implies  that  /  increases  monotomcally  when  c  >  0  Addi¬ 
tionally  /  increases  fast  enough  so  that  the  slope  of  the  secant  from  the  origin  to 


Figure  2.1.  A  re p rese n l a 1 1 ve  strain  curve  /(e)  (top)  and  the  monotone  function 
m  ( f )  which  generates  it  (liottom)  /  is  asymptotically  linear  for  small  e  At  A. 
the  slope  of  /  is  infinite  At  B,  /  is  tangent  to  the  secant  to  B  from  the  origin 
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any  given  point  on  the  curve  /(f),  is  less  than  the  slope  of  the  secant  to  the  given 
point  from  any  other  point  on  the  curve,  when  t  is  positive  for  !>oih  points 


!U\)  /U2)  /Im> 

I  •  i  )  -  -  — -  <  -  - 


-  i  1  fj 


l  |  £  •_.  >  t) 


(2  3  7  | 


-  I  -  «  1 

We  prove  this  first  for  £  ■>  >  fj,  in  which  case  m  1 1  ■  < )  >  m|>  ||,  and  thus,  /(£j)  > 
m  |i  |  )f>  We  subtrait  /(tj)  =  rn(t|)£|  from  the  latter  in  u  a  1 1  ly .  divide  by  the 
positive  i|iiantity  ?  ■_.  fj,  and  obtain  the  cfesired  result  On  the  other  hand,  when 

-  <  f  !  we  have  ill  (  *  J  )  c"  III  (  £  |  )  Then  /  (  £  .  |  <  III  ( ‘  |  ) :  a  a  II  o  /(f'.’l  >  -  fTI  ( £  j  )  £  j 

Ifv  adding  /  ( £  |  )  ==  m(£|)£|  to  the  latter  inequality,  and  dividing  by  f|  -  f2,  the 

proof  is  complete 


In  the  proof  of  (2  3  7),  there  is  the  tacit  statement  that  for  all  0  <  c  <  £0,  the 
graph  of  /(£)  lies  below  the  secant  from  the  origin  to  the  point  (e0l/(e0))  We 
state  this  separately  for  future  use  as 


/(r)  <  m  (e0 )e .  0  <  £  <  £0 


(2  3  8) 


Similarly,  we  have 


/(c)  >  m ( £ „ ) £ ,  £  >  £ 


(2  3  9) 


We  w  ill  have  occasion  to  consider  £  as  a  function  of  f 
of  /  ensure  that  it  has  a  strictly  increasing  inverse  F.  such 
inverse  has  properties  similar  to  those  of  / 

>  </)  -  /  /  *  <> 

A  I  /  )  <  1 ,  /  >  <  i 


The  above  properties 

i 


that  F ( / ) 


f 


This 


I 

t 


(2  3  Id) 


i 

i 


for  each  positive  M  <  I  there  is  a  unique  non-zero  value  of  /  which  satisfies 


r*  m 


mV 


2.4  Regularly- Vary ing  Functions. 


The  presentation  in  this  section  is,  for  the  most  part,  ident  al  to  that  given 
by  Pipkin  (2.2|. 

Real  material  response  functions  such  as  G (t )  vanish  for  f  <  0  and  are  strictly 
positive  for  t  >  0  Because  of  the  wide  variation  in  both  the  values  of  G  and  the 
time  scales  over  which  these  changes  take  place,  they  are  usually  specified  by  graphs 
on  a  doubly-logarithmic  scale  [2.3].  Let  p(f)  be  the  slope  on  a  dou  bly -logarit  h  m  ic 
plot,  for  t  >  0,  of  a  differentiable  response  function  /(f): 

p(f)  =  d(ln /(f))/d(ln  f)  =  f/'(f)//(f),  (2  4  1) 

where  /'  is  the  derivative  of  /.  We  note  that  p(t)  =  p  if  /  =  ctp  In  the  general 
case,  by  integrating  (2.4.1)  we  obtain 

ln(/(cf)//(f))  =  J  p(tx)(d  x  /  x ).  (2  4  2) 

Then  if  p(t)  approaches  a  constant  value  p,  say,  as  (  -♦  oo  ,  the  integral  in  (2  4  2) 

approaches  p1l(c),  and  thus 

f(ct)  ~  cF/(t),  l  —  oo  (2  4  3) 

H  egularly -varying  functions  are  those  functions  for  which  f{cl)/f(t)  approaches  a 
finite  non-zero  limit  as  f  —  oo  ,  and  for  any  such  function  the  limit  is  necessarily  of 
the  form  cp  as  a  function  of  the  parameter  e,  with  -  oo  <  p  <  +  oo  (see  Feller  (2  4)) 

For  such  a  function  we  write  /*(/)  =  p  and  say  that  p  is  the  power  of  / 

If  equation  (2  4  3)  holds  with  p  =  0  the  function  is  tlo  «r  Ip  ■  ns  rym  g  We 

use  I.  ( t )  (to  suggest  a  logarithm)  to  stand  for  any  otherwise  unspecified  slowly- 
varying  function  Then 

L(cf)  -  L(l).  I  —  oo  ,  (2  4  4) 
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and  any  regularly-varying  function  /(l)  can  be  expressed  in  the  form 

/(<)  =  L(t)tp  (2.4.5) 

For  many  purposes  L  can  be  treated  as  if  it  were  constant,  even  though  it  may 
diverge  to  zero  or  infinity  For  example,  powers  and  products  of  regularly-varying 
functions  are  also  regularly-varying,  with  the  obvious  exponents.  In  particular,  any 
function  which  approaches  a  non-zero  constant  value  is  slowly-varying. 

The  assumption  that  pfl)  (in  equation  (2.4.1))  approaches  a  limit  implies  not 
only  that  /  is  regularly-varying  but  also  that  /'  is  regularly-varying  (unless  /  is 
merely  slow  |y -vary  in  g  )  If  p(<)  — •  p  then  from  (2.4.1), 

f'-pf/t,  Pt  0 

1  —  oo  .  (2.4.6) 

/'  =«(//<).  p  —  0 

The  exception  for  p  =  0  occurs  because  when  /  is  approaching  a  constant,  f1  may 
be  approaching  zero  much  faster  than  1/1. 

On  the  other  hand,  when  /  is  regularly-varying  its  indefinite  integral  /j  is 
always  reg  u  la  r  ly  -  va  ry  in  g  [2.4).  When  /  j  diverges  as  1  — ►  oo  , 

[‘  tf(t) 

flit)--!  /(r)rfr~- - -  p  =/>(/)>-  1 ,  1-00,  (2.4.7) 

'  0  ( P  +  '  ) 

just  as  if  /  were  actually  a  power  When  p  <  -1  the  integral  (2  4.7)  approaches  a 
constant  and  the  tail  of  the  integral  is  regularly-varying: 

fx  1/(1 ) 

I  f  ( r )  d  t  ~ - ,  p  <  —  1  ,  1— >oo  (2  48) 

i  IP  +  1 1 

1  lie  result  analogous  to  (2  4  6b)  is 

If  H)  -  «(/,).  p  =  -  1.  l-oo  (2  4  9) 

The  presentation  above,  following  Pipkin  |2  2),  is  useful  also  for  consideration 
of  functions  whuh  are  reg  u  lar  ly  -  var  y  in  g  at  the  origin  For,  /(l)  vanes  regularly 


:« 


at  0  if  and  only  if  / (1  / 1)  varies  regularly  at  oo  (Feller  [2.4]).  Thus,  these  results  apply 
with  t  — *  oo  replaced  by  t  — *  0,  for  functions  which  vary  regularly  with  power  p  at 
the  origin.  In  particular,  the  change  of  variables  t  — ♦  1/r  in  equations  (2.4.7) 

and  (2.4.8)  returns  the  same  pair  of  integral  results  for  t  — >  0. 

We  add  here  a  result  that  is  needed  later  for  the  convolution  of  two 
regularly-varying  functions.  If  /  and  g  vanish  for  t  <  0,  and  are  regularly-varying 
as  t  —  oo  (t  — >  0)  with  powers  /•’(/)  =  p  and  P  (g)  =  q,  then 


/  f(t  ~  r)g{T)dr  ~  - - — - —  tf(t)g(t)\  p,q>-\;  <—oo(<-*0) 

•o  (P  +  fl  +  1)! 


2.5  Material  Response  Functions. 


In  the  linear  theory  of  viscoelasticity  our  stress  relaxation  function  G(t)  gives 
the  time  history  of  the  stress  (per  unit  mass  density)  in  a  viscoelastic  material 
which  is  subjected  to  a  unit  step  in  strain,  e,  at  time  t  =  0  (see  Pipkin  [2.5]). 
In  spite  of  our  normalization  with  respect  to  mass  density,  we  will  refer  to  G  as 
the  modulus.  For  the  nonlinear  theory  defined  by  equation  (2.1.5)  the  stress  is  G(t) 
when  /(«:(•,<))  =  H  (1 ) ,  the  Heaviside  unit  step  function  defined  to  be  zero  for  t  <  0 
and  one  for  t  >  0.  The  modulus  is  necessarily  zero  for  t  <  0,  since  the  current  stress 
cannot  depend  upon  the  future  strain;  the  input-output  relationship  between  the 
strain  and  the  stress  is  causal.  With  this  in  mind,  we  write  the  convolution  (2.1.5) 
as: 

a(x,t)  =  J  G(t  -  r)df  (e(x,T)).  (2.5.1) 

-  00 

We  find  it  convenient  to  make  operational  use  of  such  convolutions  and  we  use  the 
following  notation  for  convolution  with  respect  to  the  time  variable: 

U'*9)(t)  =  j  g(t-r)df(r).  (2.5.2) 

-  00 

The  limits  of  integration  extend  to  ±  oo  in  the  general  case.  We  recall  that  con¬ 
volution  is  commutative  and  associative,  and  that  it  also  commutes  with  time  dif- 
ferentation : 

(/*?)'=/'*?=/*?'■  (2-5.3) 

In  formal  operational  use  of  the  convolution,  we  use  whichever  of  these  is  convenient 
at  the  time.  With  this  notation  and  these  propeities  of  the  convolution,  we  write 
the  nonlinear  constitutive  law  in  terms  of  the  modulus,  once  and  for  all  as: 


a(x,t)  =  (G  '  *  /(e))(x,  t). 


(2.5.4) 


We  consider  moduli  G  which  have  a  finite  value  G0  at  t  =  0  and  which  decrease 


to  a  strictly  positive  equilibrium  value  Ge  as  t  — *  oo  .  In  fact,  we  assume  that  G  is 
completely  monotone  for  t  >  0;  i.e.,  G  is  positive  and  decreasing,  and  it  possesses 
derivatives  of  all  orders  which  satisfy: 

(- l)n6,(n|(0  >  0  ;  n  =  0,1,2 _ ;  t  >  0.  (2.5.5) 

Since  Ge  is  nonzero,  the  materials  we  consider  are  viscoelastic  solids  (2.5],  It  is 
convenient  to  allow  the  elastic  limit ,  in  which  case  G(t)  =  G0H(t)  -  GeH(t). 

These  properties  of  G  are  enough  to  ensure  the  existence  (2.5)  of  a 
compliance  J  (t)  which  vanishes  for  t  <  0,  has  a  com  pletely-monotonr  derivative 
J  '.(t)  for  t  >  0,  and  which  satisfies: 


{J1  *  G)(t)  =  mt)  (2.5.6) 

The  compliance  takes  a  finite  jump  to  J0  =  1  /Ga  at  t  -  0  and  increases  monoton- 
ically  thereafter  since  j'+  >  0.  For  solids  J  has  an  equilibrium  value  Je  =  1  /Ge 
which  it  approaches  as  t  — ►  oo  .  The  units  of  J  are  the  inverse  of  those  of  G\  in  our 
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case,  J  has  units  of  (velocity)  .  Whenever  we  write  specific  forms  of  G  or  J ,  we 
mean  their  values  for  t  >  0,  and  leave  implicit  that  they  vanish  for  t  <  0. 

By  operating  on  equation  (2.5.4)  with  J  \  using  (2.5.3)  and  (2.5.6),  we  obtain 
the  constitutive  law  in  terms  of  the  compliance: 

f{e(x,t))  =  (J  '*  a)(x,t),  (2.5.7) 

where  we  have  used  the  fact  that  con volut ion  with  the  Heavisidestep  function  is  just 
an  integration.  Throughout  this  work  we  will  use  both  versions  (2.5.4)  and  (2.5.7) 


of  the  constitutive  law. 


Real  material  properties  might  be  stated  in  terms  of  only  one  response  func¬ 
tion,  the  modulus  or  the  compliance.  We  now  present  some  relations  between  the 
asymptotic  behaviour  of  G  and  J,  for  the  cases  of  t  — >  0  and  t  —>  o o.  We  sometimes 
assume  that  the  difference  between  these  functions  and  their  initial  or  equilibrium 
limits  is  regularly-varying.  We  do  so  for  t  — ►  0  and  take  for  the  modulus: 

G(l)  =  G0(l  -  xi(0),  (2.5.8) 

where  Xi  is  a  positive  regularly-varying  function: 

Xl(0  =  L(t)tp  ,  0  <  p  <  1,  f  -  0.  (2.5.9) 

The  power  p  is  restricted  to  these  values  since  G  must  be  completely  monotone. 

Additionally,  x  l  must  be  completely  monotone.  Therefore,  by  equation  (2.4.6) 
(for  p  /  0),  so  must  L(f)fp_1.  Since  the  product  of  two  completely  monotone 
functions  is  again  completely  monotone  (Feller  [2.4]),  any  completely  monotone 
slowly-varying  L(t)  will  do.  One  such  function  is  |ln(<)|.  nr  ’itionally,  for  p  /  0, 

L  =  1  renders  xi#  completely  monotone.  If  p  =  0,  then  L1  must  be  completely 

monotone. 

Similarly,  we  let  the  compliance  have  the  form: 

J(t)  =  J0(l  +  X2(0).  (2.5.10) 

and  seek  to  relate  X2  to  Xi-  We  use  equation  (2.5.6)  and  the  relation  J0G0  =  1  to 
obtain  a  Volterra  integral  equation  for  X2: 

X2  =  Xi  +  Xi'*X2-  (2.5.11) 

This  is  asymptotically  satisfied  with  X2  ~  Xi-  For,  according  to  equation  (2.4.10), 
we  then  have  P(xi  '  *  X2)  =  2  (x  1 )  =  2p,  and  the  convolution  is  of  higher  order 
than  x  1  ■  So, 


J  (t)  ~  J  0(  1  +  Xi  (<)). 


t  -  0, 


(2.5.12) 


and  we  see  that  G  and  J  are  approximately  algebraic  reciprocals  for  small  t. 


For  t  — »  oo ,  we  present  some  results  obtained  by  Pipkin  [2.2].  These  were 
gotten  with  the  use  of  the  Tauberian  theorems  for  Laplace  transforms  [2.4]  to  relate 
the  time  functions  J  and  G  to  their  transforms  in  the  asymptotic  limit.  There 
are  two  situations  to  consider  which  depend  upon  the  time-dependent  apparent 
viscosity  [2.5]: 


n(0  = 


/  [ G(T)-Ge]dT . 

o 


(2.5.13) 


This  integral  exists  for  all  finite  t  for  the  moduli  we  consider. 


We  consider  first,  the  cases  in  which  p(t)  is  not  slowly-varying.  Here  we 
restrict  attention  to  moduli  for  which  G  —  Ge  is  regularly-varying  and  write 


G(t)  —  G e(  1  +  X  (0)  > 


(2.5.14) 


w  here  now , 

X(l)  =  L(t)rp,  0<p<l.  (2.5.15) 

Cases  in  which  p  >  1  are  excluded  from  this  class  since  we  are  considering  viscosities 
which  are  not  slowly-varying.  Complete  monotonicity  of  G  requires  that  p  be  non¬ 
negative.  The  result  is  again  that  J  is  approximately  the  algebraic  reciprocal  of  G: 


J  (t)  ~  Je(  1  -  x(0).  t  -+  oo 


(2.5.16) 


Now,  when  rj{t )  is  slowly-varying,  including  all  cases  for  which  the  integral  in 
(2.5.13)  converges  at  infinity,  the  response  functions  have  the  forms 

G(t)  =  ce(l  +  xi(0),  J(‘)  =  -Ml  -  X2(0).  (2.5.17) 

where  we  do  not  require  xi  and  X2  * 0  k>e  regularly-varying  .  The  relation  between  Xi 
and  x 2  is; 

/  X2iT)dr  ~  f  Xi  (r)dT  =  ~P~  ,  t  oo.  (2.5.18) 

0  0  G  g 
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This  does  not  imply  that  \2  ~  X  l  •  However,  if  t]  (  oo  )  <  oo  ,  the  relation  (2.5.18) 
becomes  an  equality  in  the  limit  (2.2],  As  a  simple  example  of  this  case,  let  G 
and  J  be  given  by  (2.5.17)  for  all  t  >  0,  and  let  Xi(0  =  exp{  —  f/r}.  With  the 
use  of  (2.5.6),  we  find  that  X2(0  =  ( 1  /2)  exp  {  —  </2T  } ,  which  is  clearly  not  asymp¬ 
totic  to  xi-  However,  both  integrate  to  T  in  the  limit,  and  the  limiting  viscosity, 
r)(oo)  —  TGC,  is  finite. 


2.6  The  Constitutive  Law:  Some  Properties 
<fc  the  Quasi-Elastic  Approximation. 


In  this  section  we  examine  the  instantaneous  and  equilibrium  elastic  behaviour 
predicted  by  our  constitutive  law  for  a  nonlinearly  viscoelastic  solid.  We  associate 
with  this  behaviour  two  wavespeeds  which  are  relevant  in  this  work.  We  then  present 
an  approximation  to  the  viscoelastic  constitutive  law  which  treats  the  material  as 
elastic  with  a  time-varying  modulus  and  compliance.  For  cases  of  a  monotone  input 
function  in  the  constitutive  law  and,  separately,  a  non-negative  input,  we  present 
simple  bounds  on  the  output  function. 

Consider  a  perturbation  strain  history  at  a  particular  location  x  which  vanishes 
for  t  <  0,  jumps  to  e0  at  t  =  0+,  approaches  ee  as  t  — ►  oo,  and  which  is  continuous, 
but  otherwise  arbitrary  for  t  >  0.  From  the  constitutive  law  (2.1.5),  the  stress  is: 


c[x,t)  =  f(e0)G(t)  +  J  G(t  -  r)fT(e(x,  r))  dr. 

■v 

For  t  <  0,  the  stress  is  zero.  Immediately  after  the  jump,  it  has  the  value: 


(2.6.1) 


o(x,0+)  =  f(e0)G0.  (2.6.2) 

On  the  other  hand,  in  the  limit  t  — »  oo  the  equilibrium  stress  from  (2.1.5)  is: 

o(r,oo  )  =  f(et)Ge.  (2.6.3) 


Equations  (2.6.2)  and  (2.6.3)  are  the  nonlinear  elastic  stress-strain  relations  which 


describe  the  instantaneous  and  equilibrium  behaviour  of  our  viscoelastic  constitutive 
law.  In  a  theory  of  wave  propagation  in  nonlinearly  elastic  materials,  the  slope  of 
a  stress-strain  curve  at  a  point  (  e ,  cr )  gives  the  square  of  the  speed  of  travel  (the 
wavespeed)  for  continuous  disturbances  having  these  values  of  stress  and  strain. 
(The  slope  of  a  secant  of  the  curve  is  the  square  of  the  speed  for  discontinuous  solu¬ 
tions  (shocks),  which  have  jumps  corresponding  to  the  values  of  the  stress  or  strain 
at  the  endpoints  of  the  secant.)  For  our  nonlinear  theory  of  viscoelasticity,  we  will 
make  use  of  the  instantaneous  and  equilibrium  wavespeeds  at  zero  strain  [1.5],  which 


I 

are  defined,  respectively,  by:  | 

i 


U]  =  lim  G 0f  '(e)  =  G0  =  J~  1 
*10 

U2  =  lim  Gef'(e)  =  Ge  =  J  ~  1 . 
*10 


< 

t 

(2-6.4)  ; 


With  the  limiting  elastic  behaviour  as  motivation,  we  now  present  an  ap¬ 
proximation  which  reduces  the  viscoelastic  problem  to  one  of  nonlinear  elasticity 
for  all  time,  with  time-varying  modulus  and  compliance.  We  call  it  the  quasi-elastic 
approximation;  it  amounts  to  approximating  the  convolution  by  multiplication  in 
either  version  of  the  constitutive  law: 


<7q(  0  =  G(t)f(eq(t))  ~  (O' *  f(eQ))(t), 
and  (2.6.5) 

/(eQ(0)  =  nt)cQ(t)  * 

where  the  z-dependence  of  the  stress  and  strain  is  left  implicit  since  it  has  nothing 
to  do  with  the  approximation.  We  call  any  such  function  eq  or  oq  which  “solves”  a 
viscoelastic  problem  based  on  the  use  of  this  approximation  a  quasi-elastic  solution. 


The  relations  (2.6.5)  become  equalities  if  the  material  is  purely  elastic.  With  this 


approximation,  we  treat  the  material  as  if  this  were  true  at  all  times  by  supposing 
that  Cl  (t)J(t)  =  1  at  each  value  of  t 


In  cases  of  a  m o n o ton ically  increasing  input  to  either  version  of  the  constitutive 
law  ,  (251)  or  (2.5.7),  we  can  obtain  sim  p le  bounds  on  the  output  function  Consider 
first  (2.5.4)  written  without  the  x -dependence  as 


(G 


r)  df ( t  (  r  ) ) 


( 2  6  G  I 


-  X 

If  £  is  increasing  m  on  o  ton  ically ,  then  so  is  /(=),  and  therefore  <  if  >  0  Now  (1 
satisfies  (le  <  (1(1  —  r)  <  Cl 0  for  -  oc  <  r  <  I,  so  we  have  the  general  bounds 


G  e  /( e  (<))<(  G  '  *  f  (e))(  t)  <  G  0f  ( ;(/)),  -  oc  <  t  <  +  oc  (267) 


Similarly,  from  (2.5.7): 


J0a(t)<(j'*o)(t)<Jeo{t),  —  oo  <  /  <  +  oo  .  (2. 68) 


If  the  input  function  vanishes  for  t  <  0,  then  the  lower  bound  on  G(t  -  r)  and  the 
upper  bound  on  J  (t  —  r)  can  be  made  tighter,  with  the  result  for  t  >  0: 


C(t)/(e(t))  <  (O'  *  /(e))  (  0  <  C0f(e(  t))  (e  =  0,t  <  0) 
J0a{l)  <  (J  '  *  <r)(  0  <  J(t)a(t)  (a  =  0,  t  <  0). 


(2.6.9) 


We  note  that  the  inequalities  for  the  bounds  in  equations  (2.6  7)  to  (2  6  9)  become 
equalities  when  the  material  is  purely  elastic  Observe  that  for  the  tighter  bounds 
given  by  equation  (2.6.9),  the  quasi-elastic  approximation  serves  as  a  lower  bound 
on  the  output  stress  and  as  an  upper  bound  on  the  output  strain.  We  also  note  that 
the  quasi-elastic  approximation  for  any  m  onoton  ically  increasing  input  always  falls 


within  th<*  general  Imumis  It  in  for  such  inputs  a  rignuroun  approximation  in  that 


bounds  mi  t  h  <■  error  (however  loosej  are  given 

U  hen  the  input  function  is  known  only  to  l>e  non  negative  the  upper  bourn 
ill  ( 2  t>  7)  ami  the  lower  liouml  hi  (2  (>  S }  are  still  valni  We  ilerive  the  latter 

(  J  '  •  a  )  I  I )  j  n  ( /  r  )  if  ./  ( r  | 
n 

J  .  CT  (  f|  a  I  a  {t  :  I ./  \  I  r  |  ,i  ? 


J.  r>\t\ 


l'2i.  1 1 


where  we  have  useii  equation  (2  2  I)  to  integrate  out  the  jump  in  ./  The  result 
obtains  since  n  is  n  on  -  n  eg  a  t  tv  e  ami  J  +  is  positive  since  it  is  completely  monotone 
The  inequality  is  strict  unless  a  |  r  |  vanishes  identically  for  all  r  <  (  Similarly,  the 
upper  bound  in  (2  C  7)  results  from  (1  '  <  0 


CHAPTER  3:  STEADY  WAVES 


3.1  Governing  Equations. 


\S  e  st  inly  steady  waves  These  are  disturbances  of  unchanging  shape  which 
propagate  w  ith  a  constant  speed  In  one  space  dimension. 

u(x,f)=u(<-x/f')  (.‘till 

describes  a  wave  of  shape  u  which  travels  with  the  constant  wavespeed  ('  VN  ith 
/'  >  0,  the  wave  travels  “to  the  right,”  in  the  positive  x-direction  for  increasing 
time  t  We  call  the  quantity  r/  defined  by: 

if  (x  1 1 )  —  t  —  x/U  (3.1.2) 


the  retarded  tim  e. 

If  the  displacement  u(x,<)  is  given  by  (3.1.1),  then  the  strain  e  and  the  veloc¬ 
ity  v  necessarily  have  the  same  form: 

«'(«»)=«  '(l)  -  Ui(ri),  (3  1 .3) 


where  the  prune  denotes  differentiation  with  respect  to  the  argument.  To  show  tha' 
this  is  also  true  for  the  stress  <r ,  we  use  equation  (2.5.1)  to  write 


o  (z ,  t) 

We  use  the  fact  that  t  —  r 
ables  r/(x,  j)  —  r  to  obtain 


f  G(t  -  r)df[i(r,(x,r)))  (3  14) 

-  OO 

i f[x,t)  —  i ;(x,r),  followed  by  the  change  of  vari- 


<7  (x,  f ) 


r  n(z.t) 


<7(i;(x,<)  -  r)  df{i(r)) 


(3  15) 


I  h.T.-f..r.>  the  stress  i.s  also  a  function  of  t h <•  retarded  tm 


In  the  rust  of  this  work  tin-  displacement  strain  and  velm  it;  .»•-  will  a*  t  h  > 
stress  will  always  In  functions  of  the  retarded  time  anil  for  n  <>  t  a  1 1  ■  n  a  I  simple 
it  v  w  ill  In-  1 1  «*n  '>  ti‘>l  li>  u  |  /  )  ft  i  T  h  e  d  er  i  va  1 1  v  es  of  l  hese  f  u  n<  t  m  n  >  w  i  II  I.  e  d  e  ri  "  t  e  d 
in  n  '  |  I )  a II  'I  So  oil 

\\  i-  n  . .  w  tie  r  iv  i*  a  n  on  I  in  <■  a  r  o  r<l  in  a  r>  i  n  t  eg  ro  -d  i  lie  run  t  ta  I  equation  a  u  <1  an  initial 
.  nii'lition  w  IikIi  describe  tliu  propagation  of  steady  waves  of  strain  in  tt  o  n  li  in- a  r  I  y 


\  IM'nr 

Iasi  1 1  solids 

w  h  n h  uc 

cupy  all  of 

one  -d  i  m  e  n  sio  n  a  1  sp  ac  i 

e  for 

all  time  A  s 

hr  si  si 

1  e  p  .  we  i]  se 

q  u  a n  I  it  les 

like  (3  1  1  ) 

in  the  momentum  anil 

1  i  o  III 

pa 1 1 b ila  y  eq  ii  a 

lions 

|2  1  3)  a  nil 

(2  1  4)  to 

relate  the 

stress  to  the  strain 

tt  ,th 

the  not  at  ion  a 

i  o n v  u n t  ion s  aliovu  w  e  olitain 

o'(t)  =  i'2e'(t)  (if  1  <») 

Wei  all  that  our  stress  and  strain  are  perturbations  on  a  quiet  state  of  the  material 
In  the  distant  past  of  physical  time,  and  at  places  far  ahead  of  the  current  position 
of  the  disturbance,  the  strain  and  stress  are  zero  In  terms  of  the  retarded  time,  this 
imposes  the  conditions 

i  i  m  e  (  ( )  =  0  and  lim  <7  (  < )  =  0  (317) 

I  —  -  00  t  —  -  cc 

I  sing  (3  1  7),  we  integrate  (3  16)  from  -  oo  to  t  and  obtain 

<7(0=  l'2e{  t)  (3  1  -8) 

We  now  substitute  for  the  stress  in  (3.1.8)  from  the  constitutive  law  (2.5  4)  The 
result  is 

e(t)  =  U~2(G'  *  /(e))  (  I)  (3  19) 
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fc,"*  'lT*  . 


If  i  ii  s  t  *■ a d  we  iim1  f»pi  ,4 1 1<  i  n  { 3  I  k  )  f.  i  r  l  h  >trrs^  in  the  *  <  >  n  s  1 1 1  u  1 1  xe  law  (2  7  i  w  * 


•  >  t  >  1  .till 

/MM)  /  -  )  (  / 1  [  :i  1  l  n  i 

I-  i(  n  .1 1  in  11  s  |  3  1  'M  »ii. I  (.'1  I  |ii  I  ar . jiin.«l'- nt  s  t  a  t  in  '■  n  t  s  "f  t)i.-  int'i'r.  <1 1  ff  •  •  r  1 1 1 1 .» I 

T'l.f  t  i.>ii  w  Ii  1 1  Ii  describes  t  Ii  <■  p  r<>  p a g  a  1 1< .  n  ■ .  f  -  t  <■  .ft  y  strain  waves  in  t  ii »-  in  a 1  ••  r  i  a  I  s 
we  •  "  11  s  til  i>  r  We  •  a  II  a  f  u  11 1  1 1"  n  •  (  f  |  a  filiation  <4  ■  1 11  r  problem  if  it  > .« t  is h ns  e<|  n  a 
( 1 '  ■  11  ( .'i  I  o|  ati'l  t  Ii  initial  •  •  •  11  •  i  1 1 1  •  >  ti  ( .')  I  7|  If  such  a  solution  exist-  the  .  ..rre 

s  j"  1 11 ')  1  ii  1’  stress  i-  "litaine'l  'Iirei  t|\  from  ei|  11  a  t  n .  n  (It  1  s  |  the  \e|...  itv  f  r  •  •  1 1 1  (3  I  3 

I'  lie  '1  |s  ||  |a>  e  III  e  II  t  Is  It  |\  e  ||  I .  \ 

ti  I  f  I  i  f  i  |  r  I  d  r  /  <  x  ( :t  I  I  1  I 

.X 

I  "i  the  f  11 11 1  t  inn  s  we  consider  (§1)2  1  and  §§2  2).  this  integral  is  w  ell-deti  tied  f"t  all 

tl  II  It  e  ( 


3.2  Steady  Wave  Solution*  in  Linear  Viscoelasticity. 


Wo  prepare  to  study  the  steady  wave  problem  for  nonlinearly  viscoelastic  ma¬ 
terials  by  first  considering  vhe  linear  problem  in  which  the  strain  curve  is  /(f)  =  e 

-'(f)-  f  *  (  J  '  *  f  )  (  / ) ,  f  ( -  oe  )  =  0 .  -  oc  <  I  <  4  oc  (3211 

l  liis  system  has  continuous  solutions 

f  ( ( )  -  er  ( 1  ~  *'  * ,  -  oc  <  f  <  4  x  (3  2  21 

where  I.,  is  an  arbitrary  time  shift  and  for  which  the  “rise  time"  l/r  satisfies 

/■'  * rJ(r )  -  1  (3  2  3  1 


w  h  ere 


rJ(r) 


t  Tl  dJ  (/) 


2.r> 


(■•f  2  4) 


1  If  integral  in  |  t  2  I  is  the  I.  a  p  Ia<  S 1 1«- 1 1  j*-s  transform  ..f  J  with  l  h  **  real  t  r  a  u  sfo  r  in 

•  i  •  i  i  i  I  *  1  i  '(  n  i  \  t  n  1 I  v  .1  is  t  ii  ••  I  a  |>  I  a i  •-  transform  <  >  (  J  ,  1 1  i s  w  «•  1 1  il  ••  ti  n  •- . I  a  1 1  |j 

f  •  ii  ■  i  f. .  i  lii-  ■  f  in  (. !  i  a  ii  '  «-s  ./  fa  .  i  i,  ii  s  i  <  I  c  r  I  'I  i, 2  ‘F  i 

W  •  I.  •-«  show  I  II  a  l  III  I'll'  |s  a  111111(111-  I  <•  lat  If  f  II  sli  I  [I  between  all  pos.it  IV  «-  value' 

I  '  i  n  I  tii*  j  ■  >  •  s  1 1  i  \  *■  v  a  In  rs  . .  f  /  >*  iii<  i>  satisfy  equation  |  .1  2  3  )  I  s  i  n  g  the  st  a  mi  a  r  ii 

i  •  -.it  t  •  :  lli*  1 1  in  1 1  s  ■  (  s  ii  i  (j  t  r  a  li  si'  >  I  Hi  s  w  .  have 

i  i  in  r  ./  I  r  )  1 1  in  ./  |  /  I  ./  I 

*  •  a  r  •  n 

.'I  2  '<  ! 

1 1  Hi  r  .1  |  r  |  llin  ./  |  /  |  I  , 


h  Iff.-  vs.  lias.-  us.-.)  equations  I  2  (>  i)  l<>  r  <■  I  a  l  «•  tlf -si-  limits  t  <  >  tin*  i  n  st  a  n  t  a  neo  u  s  and 

i"i  >i  ill  l>  r  in  in  v*  as  *sj . i  s  \i,»  ,{  J  .>  (i  implies  that  r./(r|  is  (completely)  rn  on  <.  to  i  e 

a  i  I  n  m  l  i'  f  u  of  r  2  I  ,  thus  t  |i  .■  relationship  between  r  g  (0.  oc  )  and  f '  g  |  f, .  V  0  )  is 
•  if  I-  mu  <'  ii  nt  iii  umis  solutions  to  (3  2  I  |  of  tin*  form  (3  2  2)  therefore  exist  only 

f'  i  ■  h  values  of  I 

W  e  now  consider  discontinuous  solutions  to  equation  (3.2  1)  We  call  such  a 
solution  a  ehork  if  it  possesses  a  finite  jump  discontinuity  of  the  form  (2  2  2)  If  the 
solution  is  continuous  but  its  derivative  has  a  jump  discontinuity,  then  we  call  it  an 
acceleration  Mine,  if  such  a  discontinuity  appears  in  a  higher  order  derivative,  the 
solution  is  a  higher  order  acceleration  uane  It  is  well-known  (3  1]  that  shocks  can 
f"  r  in  i  n  h  n  1 1  e  time  in  elast  ic  m  a  te  r  i  a  Is  w  hi'  It  obey  a  nonlinear  stress-strain  relation 
>  i  m  i  la  r  Iv  the  ei|  u  a  I  ion  s  governing  linearly  e  last  ir  media  ran  support  traveling  shock 
solutions  provided  the  discontinuity  is  in  the  initial  data  of  the  problem  these 
disturbances  travel  at  the  constant  wavespeed  f  where  a  —  (  ~f  is  the  linear 
si  re  ss-st  rain  relationship  We  now  show  that  equation  (3  2  I)  governing  steady 
waves  for  the  theory  of  linear  v  iscoelasticitv  admits  non-negative  discontinuous 
solutions  among  the  class  of  functions  possessed  of  a  Laplace  transform  if  and  only 


.1 


if  tin-  material  is  purely  elastic  Suniciemy  is  immediate,  since  the  problem  reduces 
t"  the  linear  elastic  problem  Necessity  is  obtained  m  the  following  We  seek  a  non- 
negative  solution  -  which  vanishes  for  f  <  (l  in  the  class  of  functions  considered 
"  e  t  ran  sform  eq  ua  t  mu  |  3  2  1  )  using  equation  (3  2  d  )  a  Ion  g  with  the  fact  that  the 
transform  of  the  convolution  is  the  product  of  the  transforms  of  the  convolved 
functions  ami  » I  >  t  a  in 


|  .1  )  =  /  ">  ./(  s 


(■'I  2  (,) 


w  here  s  is  t  h  e  complex  t  ra  n  sform  variable  for  r  >  0  and  not  identically  zero,  e  (  s  )  is 
a  n  analytic  f  u  net  i  o  n  of  s  which  is  real  and  posi  1 1  ve  for  a  II  re  al  n  >  r0.  for  some  r0  >  0 
Therefore,  n ./  |  a )  -  I'  '  on  the  segment  s  >  r(,  Its  unique  analytic  continuation  to 
the  whole  complex  plane  is  *  J  ( * )  =  I  ~  ~  .  this  inverts  to  give  7(f)  =  U  ~  *  H  (< )  .  We 
conclude  that  non  negative  steady  shocks  and  acceleration  waves  which  vanish  on 
the  interval  (  -  oc  ,  0)  and  possess  a  Laplace  transform  exist  in  the  linear  theory  of 
viscoelasticity  only  for  purely  elastic  materials 


3.3  General  Characteristics  of  Nonlinear  Steady  Waves. 


We  now  consider  non-negative  solutions  to  equation  (3  19)  governing  non- 
linearly  viscoelastic  steady  waves  subject  to  the  initial  condition  (3  17)  In  this 
-,ci  tion  we  assume  the  existence  of  these  solutions  As  motivation  for  our  investi¬ 
gation  w  <■  lirst  take  a  geometrical  approach  to  this  system  W e  then  derive  a  set 
of  necessary  conditions  on  a  solutions'  behaviour  These  conditions  are  presented 
in  a  list  of  general  •  It  a  r  ac  ter  is  t  ic  s  of  lion  lin  early  viscoelastic  steady  waves  We  call 
them  firnpr  rtir.i  In  some  cases,  they  an1  statements  of  nonexistence.  The  proofs  of 
these  properties  follow  their  summary  We  then  use  these  results  along  with  some 
results  from  the  previous  section  to  estimate  the  graphs  of  solutions 


JiVLV.V  V  A"  AW. 


Consider  equation  (3  1.9)  written  in  the  form 

U2e  =  G  '  *  f(e).  (3.3.1 ) 

From  equation  (2.6.7),  we  have  that  for  m  on  o  to  n  ically  increasing  strains  e,  the 
stress  '  *  /(f)  is  bounded  below  by  Cef(e)  and  above  by  G0f(e).  On  a  graph 
of  stress  vs  strain,  this  stress  always  lies  between  its  equilibrium  and  instanta¬ 
neous  curves  On  such  a  graph,  the  left  side  of  (3.3.1)  is  a  line  with  slope  U~ 
passing  through  the  origin;  the  wavespeed  V  is  a  parameter  in  the  problem.  Con¬ 
sider  Figure  3.1  and  recall  the  definitions  of  the  equilibrium  and  instantaneous 
w avespeeds,  {',  and  U0.  given  in  equation  (2.6.4).  We  will  see  that  the  value  of  U 
in  relation  to  and  V 0  determines  the  nature  of  the  solution  e.  From  the  graph 
in  Figure  3.1  and  equation  (2.6.3),  we  expect  that  a  solution  has  an  equilibrium 
value  where  the  line  intersects  the  lower  curve.  Similarly,  with  equation  (2.6.2),  a 
solution  should  exhibit  a  non-zero  instantaneous  (discontinuous)  response  if  the  line 
intersects  the  upper  curve  other  than  at  the  origin.  These  properties  are  among 
those  of  our  desired  solutions.  We  prove  that; 

1  If  f(  /)  is  a  steady  wave  solution  such  that  c(  t)  — *  ee  >  0  as  t  — ►  oo  ,  then  the 
equilibrium  value ,  ee,  is  the  unique  non-zero  solution  to 

/(ef)  =  U2Jtee.  (3.3.2) 

If  (’  >  Uf.  then  €  e  is  strictly  increasing  as  a  function  of  U  . 

2  There  are  no  non-negative  solutions  possessed  of  a  non-zero  equilibrium  value 
if  V  <  V , ,  hence,  the  only  non-negative  monotone  solution  at  these  wavespeeds 
is  the  trivial  solution  e(f)  =  0.  If  V  >  Ue,  then  a  non-trivial  non-negative 
monotone  solution  has  a  finite  equlibrium  value  e,  >  0. 

3  For  all  values  of  U  ,  if  c  ( 7^ )  =  0  for  some  finite  T,  then  e  =  0  for  all  t  <  T . 


4.  Any  jump  discontinuity  satisfies  the  shock  condition: 


J  o 

~=U2J0,  (3.3.3) 

A  e 

which,  for  our  strain  curves  /,  means  that  any  steady  shock  travels  at  a  wave- 
speed  U  >  U0. 

5.  If  U  >  U0,  then  any  non-trivial  non-negative  solution  e  cannot  be  globally 
continuous.  It  must  exhibit  one  and  only  one  jump  discontinuity.  This 
jump,  at  t  =  0  (say),  is  from  zero  to  s0,  the  unique  solution  to 

/{£„)  =  U2J0e0.  (3.3.4) 

For  all  t  <  0,  e(  t)  is  identically  zero.  As  a  function  of  U ,  e0  is  strictly  in¬ 
creasing.  Additionally,  e0  <  £e,  with  equality  if  and  only  if  the  material  is 
elastic. 

6.  If  e(  <)  is  a  non-trivial  continuous  solution  which  vanishes  for,  say,  /  <  0, 
then  U  =  U0  (i.e.,  U0  is  the  wavespeed  for  acceleration  waves). 

We  first  prove  Property  1.  It  is  convenient  to  use  the  equivalent  formulation 
of  the  problem  given  in  (3.1.10)  written  as: 

/oo 

e(f  -  t)  dJ(r).  (3.3.5) 

o 

If  e(t)  — *  £e  >  0  as  t  — ♦  oo,  then  we  obtain  (3.3.2)  as  the  limiting  form  of  (3.3.5), 
since  /0°°  dJ  (r)  =  Je,  according  to  (2.2.5).  From  the  properties  of  /  (§§2.3),  non- 
zero  solutions  to  (3.3.2)  are  unique  and  increase  strictly  with  U  when  U  Je  >  1,  i.e., 
when  (I  >  Ue.  D 

We  remark  that  one  could  also  pose  our  problem  as  a  search  for  steady  waves 

having  a  given  non-zero  equilibrium  value,  since  for  a  given  a  compliance  J  and  a 

particular  ee  >  0,  the  squared  wavespeed  U2  is  uniquely  determined  from  (3.3.2). 
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Figure  3.1.  The  instantaneous  and  equilibrium  stress-strain  curves  and  the 

line  U2e. 

To  prove  Property  2,  we  observe  that  the  unique  equilibrium  strain  from  equa¬ 
tion  (3.3.2)  is  zero  when  U  <  Ue  for  our  strain  curves  /.  Thus,  there  are  no  non¬ 
negative  solutions  with  non-zero  ee  at  these  wavespeeds.  Since  a  solution  e  vanishes 
as  t  — »  —co,  it  can  be  monotone  and  vanish  as  I  — »  -t-  oo  if  and  only  if  it  is  identi¬ 
cally  zero.  Now.  suppose  that  e(l)  is  a  non-trivial  non-negative  monotone  solution 


to  (3.3.5)  for  some  value  of  U  >  Ue.  We  use  the  upper  bound  on  the  convolution 
given  in  (2.6.8)  to  write 

f{s(t))  <  U2Jee(t).  (3.3.6) 


Recall  (§§2.3)  that  m  (e)  =  f  (z) / £  is  strictly  increasing  from  m  (0)  =  1.  As  such,  it 

has  a  strictly  increasing  inverse  m_1,  defined  for  m  >  1,  such  that  e  —  m_1(m). 

with  m_1(l)  =  0.  We  use  this  in  equation  (3.3.6)  to  obtain 

e(<)  <  m  l(U2Je)  for  all  t.  (3.3.7) 

If  e(t)  is  a  non-trivial  non-negative  monotone  solution,  equations  (3.1.7)  and  (3.3.7) 
imply  the  existence  of  a  finite  limit  se  >  0  for  t  —>  oo  .  □ 

Property  3  follows  by  considering  e(T)  =  0  for  some  finite  T.  The  integral 

in  (3.3.5)  must  vanish;  with  dJ  >  0  and  for  e  non-negative,  we  conclude  that  e(t) 
must  vanish  for  all  t  <  T.  □ 


We  now  prove  Property  4.  Suppose  e  has  a  countable  number  ofjump  disconti- 
nuities  Ar;  =  e(f(+)  —  e(l.  ),  for  i  =  1,2,  ...,  where  the  ti  are  labeled  in  increasing 
order,  and  the  superscripts  +  and  —  denote  limits  to  the  point  from  above  and 
below,  respectively.  Using  equation  (2.5.3),  we  write  (3.3.5)  as: 


f{e(t))  =  U2f  J  (t  —  t)  de(r). 


(3.3.8) 


Consider  the  fc-th  jump;  for  t  <  tk,  we  have  from  (3.3.8): 

/t 

J(t-  T)e'(r)dT,  t<tk, 


(3.3.9) 


i  <  k 


for  which  we  recall  that  J  is  zero  when  its  argument  is  negative.  Similarly, 

/(«•(<))  =  U2Y,  J(t  -  ti)&ei+  U2  j  J(t  -  T)e'(T)  dr,  t  k  <  t  <  t  k+ l  {3 .3  AO) 


i<  k 


In  the  limit  as  t  |  t k  in  (3.3.9)  we  have  that 

f(e(rk))  =  U2Y,J(t;  -  e,)Ael  +  U2  f  *  J(1A:  -  r)e  '(t)  dr, (3. 3.1  1 ) 


i<  k 


and  as  t  (  tk  in  (3.3.10): 

f  ** 

!(e[l+))  =  U2J0  Aek  +  U2Y,J(t+  -  t,-)  Ae,  +  U2  J  (t+  -  r)e'(r)dr,(  3.3.12) 


where  we  have  used  J0  =  J(0+).  We  subtract  (3.3.11)  from  (3.3.12),  and  use  the 
continuity  of  J  for  t  >  0  to  cancel  the  summations.  The  integrals  also  cancel,  since 
the  integrands  differ  only  at  the  single  point  t k.  We  are  left  with: 
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U2J0  = 

Aek 


(3.3.13) 


which  is  the  desired  shock  condition  for  the  £-th  jump  discontinuity.  Geometrically, 
Af/Ae  is  the  slope  of  the  secant  on  the  curve  /(e)  connecting  the  values  of  e 
which  comprise  the  jump.  Now  when  both  e  +  and  £~  are  positive,  we  have  from 
equation  (2.3.7)  that  this  secant’s  slope  is  larger  than  one,  since  m  >  1.  If  one 
of  or  e~  is  zero  (say  e“),  then  Af/Ae  is  the  slope  of  the  secant  from  the  origin 
to  the  point  (e+,/(e+)),  and  this  slope  is  also  greater  than  one  for  non-zero  e+ . 
Therefore,  U2J0  is  strictly  larger  than  one.  This  is  equivalent  to  U  >  U0  which  was 
to  be  proven  for  any  shock  wave.  □ 

A 

To  prove  Property  5,  we  first  show  that  if  e  has  non-zero  values  for  U  J0  >  1, 
then  there  is  a  “forbidden  set”  of  values  which  it  cannot  have.  We  fix  the  current 
time  t  and  assume  there  is  a  non-negative  function  e  which  satisfies  equation  (3.3.5) 
and  whose  current  value  e(  t)  is  positive.  We  integrate  out  the  jump  in  J  and  write: 

/oo 

e(t  -  T)j'+(T)dr.  (3.3.14) 


Since  the  integral  in  (3.3.14)  is  non-negative,  we  have 


-  U2J0e(t)  >  0. 


(3.3.15) 


For  U2J0  >  1  and  for  our  functions  /,  relation  (3.3.15)  is  satisfied  as  an  equality 
by  a  unique  value  e0,  according  to  equation  (2.3.3).  As  a  function  of  U,  e0  is 
increasing  due  to  the  the  characterization  of  /  given  in  (2.3.4).  From  (2.3.9),  the 
left  side  of  (3.3.15)  is  positive  for  e  >  e0.  Values  of  e  e  ( 0,e0 )  are  forbidden  since 


WWW 


they  render  the  left  side  of  (3.3.15)  negative,  in  view  of  (2.3.8).  Therefore,  at  any 
time  Z,  either  e(Z)  =  0  or  e(  Z)  >  e0  in  (3.3.14),  when  U  2  J  0  >  1.  We  conclude 
that  the  solution  is  not  globally  continuous  (unless  it  is  identically  zero);  it  must 
jump  from  zero  in  order  that  it  have  non-zero  values.  By  Property  4,  this  jump 
must  be  to  the  value  s0.  If  e(t)  =  e0  in  equation  (3.3.14)  when  Z  =  0+,  then  the 
integral  must  vanish  at  0+.  For  non-negative  solutions  e,  we  must  have  e(r)  =  0 
for  all  r  <  0.  Furthermore,  this  jump  discontinuity  is  the  only  one  possible,  since 
equation  (2.3.7)  implies  that  a  jump  originating  from  and  ending  at  vali  es  of  e  not 
in  the  forbidden  set  must  travel  at  a  wavespeed  greater  than  the  value  of  U  which 
produced  the  jump  to  e0.  Such  another  value  for  the  wavespeed  would  contradict 
Property  4;  any  other  such  jump  cannot  belong  to  the  same  steady  wave.  Finally,  it 
is  a  trivial  consequence  of  Property  1  and  the  characterization  of  /  given  in  (2.3.4) 
that  c q  <  £ e -  If  the  material  is  elastic,  then  J0  =  Je  renders  e0  and  ee  identical  and 
this  is  not  possible  if  the  material  is  not  elastic.  □ 

To  prove  Property  6,  we  use  the  assumed  properties  of  e  to  choose  a  sequence 
of  times  {Zfc}£°=1  with  tk  j  0  such  that,  for  each  k : 


e(tk)  =  max  e(r). 


From  equation  (3.1.10),  we  have: 


/(e(Z*))  =  U2  (J  '  *  e)  (tk) . 


Since  e  vanishes  for  Z  <  0,  equations  (2.6.9)  and  (2.3.4)  imply  that 


m(e(tk))  <  U2J(tk ), 


(3.3.16) 


(3.3.17) 


(3.3.18) 


ith  equality  only  if  the  material  is  purely  elastic.  In  the  limit  as  tk  J  0  (and 


therefore  as  e  )  0),  we  have 


I  =  m(0)  <  U2J0  =  ( U/U0 )2, 


(3.3.19) 


\  > 


with  the  use  of  equation  (3.2.5).  Since  e  is  assumed  continuous,  Property  5  implies 
that  U  is  not  greater  than  U 0;  hence,  U  =  U 0.  □ 

We  now  use  the  above  properties  to  estimate  the  graphs  of  candidate  steady 
wave  solutions.  For  U  >  U0,  a  candidate  solution  is  a  shock  which  has  a  single 
jump  at  t  —  0  and  continuously  approaches  the  appropriate  equilibrium  value. 
At  U  =  U0,  we  consider  a  limiting  case  of  shock  solutions  as  the  parametric 

wavespeed  is  decreased  to  U0  to  obtain  a  candidate  acceleration  wave  (or  a  higher 
order  acceleration  wave)  which  vanishes  for  t  <  0.  For  U e  <  U  <  U0  a  candidate 
solution  is  continuous  and  increases  from  zero  at  negative  infinity  to  ee  as  t  — *  oo  . 
Since  e  must  be  vanishingly  small  as  I  -t  —  oo  ,  we  expect  it  to  be  asymptotic  to 
a  solution  of  the  linearized  equation  obtained  by  setting  f  (e)  s  e.  From  the  pre¬ 
vious  section,  we  know  that  there  arc  exponential  solutions  for  these  values  of  U . 
Thus,  our  candidate  solution  for  the  nonlinear  problem  is  asymptotically  exponen¬ 
tial  for  t  -*  -oo.  The  graphs  of  these  functions  for  all  values  of  U  are  shown  in 
figure  3.2. 


Is 


3.4  The  Q  u  a  si- E  la  at  ic  Solution. 


In  sei'tion  2  li  w<*  presented  the  q  11  asi-elast ic  approximation  to  the  constitutive 
law  We  use  it  herein  to  obtain  the  q  u  asi-e  las  t  ic  solution  Eq  for  nonlinearly  visco¬ 
elastic  steady  strain  waves  After  we  examine  its  ch  ararterist  ics,  we  discuss  its 
usefulness  as  an  approximation  to  steady  wave  solutions 

In  terms  of  the  compliance,  Eq  satisfies  the  following  algebraic  equation  at 

each  t 

f(eQ(t))  =  V2J(t)eQ(t).  (3.4.1 ) 

We  prove  these  properties  of  Eq  : 


2 


o 

1.  For  U^Jil)  >  1,  there  is  a  unique  non-zero  quasi-elastic  solution  Eq 

If  U2J(t)  <  1,  then  eQ(t)  =  0. 

2.  The  quasi-elastic  solution  has  the  correct  equilibrium  value  and,  for  U  >  Ua, 
the  correct  instantaneous  value. 

3.  Eq  is  piecewise  continuously  differentiable  and,  where  it  is  non-zero,  it  is 
strictly  increasing. 

4.  The  quasi-elastic  solution  is  an  upper  bound  on  monotone  solutions  to  equa¬ 
tion  (3.1.10)  for  nonlinearly  viscoelastic  steady  waves  which  vanish  for 
t  <  0  (say). 

5.  If  Eq  0  and  Eq  ^  are  the  quasi -elastic  solutions  corresponding  to  wavespeeds  U2 
and  Ux  for  which  (J2  >  V j,  then  cq2(<)  >  Cq^I)  with  equality  only  at  those 
tim  cs  w  here  eQ  2(  0  -£<?,(  0  -  0. 


i 


i 


4 


a 


1 


The  existence  and  uniqueness  of  Eq  (Property  1)  follow  from  the  properties 
of  our  strain  curves  /.  □ 
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For  Property  2,  the  limiting  rases  of  (3.4.1)  produce  the  results  in  equa¬ 
tions  (3.3  2)  and  (3.3.4)  for  viscoelastic  steady  waves.  □ 

To  prove  Property  3 ,  we  recall  from  sec  t  ion  2.3  that  m  (e )  =  f  (e)  /  e  is  piecew  ise 
continuously  differentiable.  Its  inverse,  in  *(hi),  defined  for  rn  >  1  is  also  piecewise 
continuously  differentiable.  Thus,  from  equation  (3.4.1)  we  may  express  Cq  direct!) 

as ; 

eQ(t)  =  m~l(U2J(t))  for  U2  J  (l)  >  1.  (3.4  2) 

The  differentiability  and  monotonicity  properties  of  m~  1  and  J  imply  these  prop¬ 
erties  of  eq .  C 

To  prove  Property  4,  we  recall  from  equation  (2.6.9)  that  the  quasi-elastic 
approximation  is  an  upper  bound  on  the  constitutive  law  for  monotone  input  func¬ 
tions  which  vanish  for  t  <  0.  Thus,  from  equation  (3.1.10),  we  have  for  a  monotone 
viscoelastic  solution  e  which  vanishes  for  t  <  0: 

f{e(t))  <  U2J(t)e(t).  (3.4.3) 

However,  equation  (3.4.1)  implies  that  U2J(t)  =  m(eQ(t)),  in  the  notation  of 
section  2.3,  and  thus, 

f(e(t))  <  m(eQ(t))e(t).  (3.4.4) 

Using  equations  (2.3.8)  and  (2.3.9),  we  conclude  that  e  <  Eq  .  □ 

The  proof  of  Property  5  follows  from  equation  (3.4.2)  and  the  monotonicity 

of  m"1.  G 

We  note  that  thegraphsofeg.determined  by  m-1  and  J,look  much  like  those 
of  our  candidate  steady  wave  solutions  shown  in  figure  3.2.  If  U  >  U0,  then  zq 
is  a  shock  which  has  a  jump  at  t  =  0,  from  zero  to  m~ 1  (U2  J  0) ,  from  which  it 
increases  monotonically  to  ee.  If  U  —  U0,  then  £q  is  a  continuous  acceleration 


wave  (of  possibly  higher  order)  which  departs  from  zero  at  /  =  0  and  increases 


monotonically  to  its  equilibrium  value.  For  the  cases  U e  <  U  <  U0,  eq  is  a  shifted 
version  of  the  above  acceleration  wave;  it  is  zero  until  some  time  l0  >  0  which 
satisfies  U"J(t0)  =  1.  It  is  obvious  that  for  this  latter  range  of  IJ  ,  such  a  quasi¬ 
elastic  solution  cannot  be  a  global  bound  for  the  viscoelastic  solution  which  never 
vanishes  on  an  interval. 

The  quasi-elastic  solution  is  attractive  for  a  number  of  reasons.  As  it  is  the 
solution  to  an  algebraic  relation,  it  is  easily  computed.  For  special  forms  of  the 
strain  curve  /,  £q  can  be  expressed  exactly  in  terms  of  parameters  of  /  and  the 
compliance  J  .  If  /  and  J  are  given  as  interpolators  of  data,  Eq  (  t)  can  be  deter  mined 
at  any  particular  time  t  from  a  numerical  nonlinear  equation  solver. 

The  quasi-elastic  solution  provides  an  upper  bound  on  the  stress  in  a  material 
experiencing  monotone  steady  wave  motion;  for,  the  stress  is  proportional  to  the 
strain  by  equation  (3.1.8),  and  eq  is  an  upper  bound  on  the  strain  (Property  4) 
In  the  next  chapter,  we  show  how  a  sequence  of  increasingly  tighter  upper  bounds 
on  the  solution  e  can  be  obtained  from  Eq.  Moreover,  for  some  shock  solutions 
(U  >  U0),  the  quasi-elastic  solution  can  be  a  rather  good  approximation  in  itself. 
We  have  shown  in  this  section  (Property  2),  that  Eq  has  the  correct  limiting  values 
at  t  =  0  and  for  <  — *  oo  .  In  Chapter  5,  we  show  further  that,  for  shock  solutions,  eq 
is  asymptotically  correct,  as  t  [  0. 
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CHAPTER  4:  EXISTENCE  Sc  UNIQUENESS  THEOREMS 


4.1  Prelude. 

In  this  chapter,  we  prove  the  existence  of  the  kinds  of  steady  wave  solutions 
we  chose  as  candidates  in  Chapter  3  and  we  examine  the  uniqueness  of  such  solu¬ 
tions  among  certain  classes  of  functions.  The  proofs  we  present  for  the  existence 
theorems  are  constructive.  We  can,  therefore,  use  the  method  of  proof  to  obtain 
a  numerical  approximation  to  the  solution.  We  build  solutions  from  m  onoton  ically 
convergent  sequences  of  functions  which  result  from  an  iteration  scheme.  Although 
monotone  convergence  is  generally  slow,  it  has  the  nice  property  that  each  member 
of  a  convergent  non-increasing  (non-decreasing)  sequence  of  functions  is  an  upper 
bound  (lower  bound)  on  the  limit  function.  We  make  these  notions  of  monotone 
convergence  precise  in  the  next  section.  In  the  following  sections,  we  present  the 
iteration  scheme  and  its  properties,  and  we  prove  the  main  theorems  on  existence 
and  uniqueness  of  steady  wave  solutions. 

4.2  Monotone  Convergence  Sc  Bounding  Sequences. 

We  denote  by  { /„ } ^°_  ,  a  sequence  of  real-valued  functions  defined  on  some 
subset  of  the  real  line.  Consider  such  a  sequence  with  members  which  satisfy: 

/l  (  0  <  f2{  0  <  <  In  (  I )  <  ...  <  /{/(<),  (4.2.1) 

at  each  t,  for  some  finite  real-valued  function  fy  (subscript  U  for  “upper  bound”), 
which  is  also  defined  where  the  fn  are  defined.  At  a  fixed  value  of  t,  the  se¬ 
quence  {  }  n  (  0  }  i  >s  a  non-dec  re  asing  sequence  of  numbers  bounded  above  by  the 
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Such  a  s<m)  ii  i'  n  i  e  ion  verges  to  a  iitimfcr  we  deiiolt-  f\  /  I  / 1  \  -  / 


varies,  a  real-valued  function  /  is  defined  on  the  set  where  the  /„  .«(.•  define. t  I  hi- 
limil  function  is  the  point  vise  limit  of  the  sequence  |  /„  }  ^  (  it  is  founded  fe|ow 
liy  each  meinfier  of  this  seeiien  c  e  and  above  by  /(  for  all  values  of  t  wher-  it  i- 
deli  ned 

/ 1  <  /_■  <  <  f  <  fr  I  I  -  I 

Naturally,  an  analogous  result  holds  for  non  -  inc  reasm  g  sequences  founded  tef.u 
by  a  lower  found  function  f/ 


■4.3  The  Iteration  Scheme. 


A  standard  technique  for  the  solution  of  integral  equations  (4  l]  is  iteration 
The  equation  to  lie  solved  is  arranged  as  an  input-output  system  into  which  an 
initial  guess  for  the  solution  is  fed.  The  resulting  output  function  is  then  used 
as  a  new  input  to  the  system  and  this  process  is  repeated,  generating  a  sequence 
of  functions.  If  this  sequence  converges  to  a  function  which  satisfies  the  original 
equation  and  any  associated  side  conditions,  the  limit  function  is  then  a  solution 
of  the  problem.  We  apply  this  technique  to  construct  a  solution  for  the  system  of 
equations  (3.1.10)  and  (3.1.7)  governing  the  propagation  of  nonlinearly  viscoelastic 
steady  waves.  We  recall  from  section  2.3,  that  the  strain  curve  /  has  an  inverse  {■ 
We  use  this  to  write  (3.1.10)  as 


I  In-  ll  i- 1  a  l  I*  •  n  st  II  •'  in  o  i- 

•  „ ,  I !  n  I  ( I'  'I  •/  '  *  -  „  H  'I )  "  1  ( l  :i  :n 

U  c  ji  r '  1 1  >■  i'll  w  i  ( (i  i  li  c  proof*  of  a  *<•  r  u*s  of  )•*  hi  in  a  *  t <■  la  t  oil  to  tin*  1 1  <*  ra  1 10  n  sell  »*iii  e 
w  li  o  h  w  ill  If  it  ti  \i  s  to  tlii*  in  a  i  n  t  li  ft  i  r<*  in  s  in  tlii-  s  u  c  < ■i-i*il  in  g  m'c  turns  Throughout,  wo 
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Proof  \\  f  hast- 

(./  '  •  f)t  f|  -  /  '  I  </■/(■  I  (4  3  -1) 

n 

Sniff  il ./  ,*  Ii,  t  li  o  in  oiiot  on  if  it  y  of  ■  implies  that  of  the  integral  □ 

Lemma  4.3.2.  If  ■  | ; )  >  -* ,  (  r  )  for  a  II  r  <  / .  with  c  ,  not  idem  ic  ally  equal  to  e  j  on 

.i  sot  of  positive  niM sure,  then  ( ./  '  *  ■:■_,)(/)  >  (J  1  *  :()(/). 

P  roof  . 

-  (J'  •  -:,)(<)  =  (./'*  (f_.  -  m))(<) 

fx 

=  J  l  )  (t  ~  r)  dJ  (t) 

o 

>  <). 

since  il .1  >  <1  and  by  the  hypotheses,  the  difference  e-<  -  £]  contributes  at  least  a 
finite  positive  amount  to  the  integral  (Z 

Now  for  the  in  p  n  t -o  u  t  p  u  t  system,  we  have  the  following  lemmas 

Lemma  4.3.3.  (a)  If  in  equation  (1.3  2)  fjn  is  monotone,  then  £oul  is  monotone, 
(h)  If  t  jn  is  continuous  except  for  a  finite  jump  discontinuity  at  l  =  0,  then  £out  ^as 


a  Unite  jump  discontinuity  at  I  =  0  and  is  otherwise  continuous 


T 


I’  roof  (  a )  I  he  m  o  no  ton  ir  ity  of  c01)l  follow  s  d  irertly  from  the  in  onoton  icitv  of  hot  h  F 
and  J  '  *  fjn  Statement  (b)  follows  from  the  continuity  of  F,  the  continuity  of  J 
{nr  I  >  0,  and  equation  (2.2.4)  for  the  Stieltjes  integral.  □ 

Lemma  4.3.4.  If  ea  and  eb  are  not  identical  and  0  <  ea  ( 7  )  <  £>J  (r)  for  all 
in  in  —  in'  —  in  ' 

r  <  /.  then  <„.(*)  <  (  t). 


Proof:  From  the  hypotheses  and  Lemma  4.3.2,  we  base 


{J1  *  =-"„)(  0  <  (■''*4,)io. 


Since  F  is  strictly  increasing, 


1 1 3.6) 


and  the  result  follows  from  (4.3.2)  and  (4.3.6).  □ 

With  the  next  pair  of  lent  mas,  we  exam  in  e  the  output  function  obtained  when 
the  input  is  related  to  the  quasi-elastic  solution  presented  in  section  3.4  The  first 
provides  sufficient  conditions  under  which  the  output  function  is  strictly  less  than 
the  input  function. 

Lemma  4.3.5.  If  Sjn  is  monotone,  vanishes  identically  for  l  <  0,  and  satisfies 
s  i  n  ( r )  >  eq{t)  for  all  t  <  t,  then  eoul(t)  <  e  in  (  f ) ,  for  t  >  0,  when  the  material  is 
not  purely  elastic. 


Proof:  By  equation  (2.3.9),  we  have  for  sin  >  Eq  : 


f{£ in)  >  w  (eq  )ein  . 


We  recall  from  (3.4.1),  that  m  (eg  (  t ))  =  U  ~  J  (  t),  so  we  h  ave, 

/(ei„(  0)  >  U2J(t)sin(t) 

>  U2(j'*ein)(t), 


(4.3.7) 


(4.3.8) 


where  the  strict  inequality  is  obtained  from  the  bound  in  (2.6.9)  for  materials  which 
are  not  purely  elastic.  Using  F ,  the  strictly  increasing  inverse  of  /,  and  the  definition 


of  s0u 1 1  wc  obtain: 

<Mn(0  >  F  (u2{J  '  *  £in)(  0)  =  ^out  (0.  (-1.3.9) 

which  is  the  stated  result.  □ 

The  next  lemma  gives  sufficient  conditions  for  the  quasi-elastic  solution  to  be 
a  strict  upper  bound  on  the  output  function. 

Lemma  4.3.6.  If  e  in  vanishes  for  t  <  0,  is  monotone,  and  satisfies  £jn(r)  <  £g(r) 
for  all  t  <  t,  then  eout(  t)  <  eq  (  /),  for  t  >  0,  when  the  m  aterial  is  not  purely  elastic. 

Proof:  By  Lemma  4.3.2  and  the  bound  in  (2.6.9),  we  have  the  following  inequalities: 

(J  '  *  cin)(  0  <  {J  '  *  £q)(  0  <  J  (  0^(3  (  0-  (4.3.10) 

So, 

*<,„»(<)  =  „)(/))  <  F(u2J{t)eQ{t))  =  eQ(t)  (4.3.11) 

by  the  monotonicity  of  F  and  the  definition  of  Eq  in  equation  (3.4.1).  □ 

Analogous  to  the  results  of  the  two  previous  lemmas,  we  can  provide  suffi¬ 
cient  conditions  for  which  the  output  function  is  strictly  greater  than  the  input 
function  and,  separately,  conditions  for  which  we  have  a  strict  lower  bound  on  the 
output.  We  consider  input  functions  which  are  related  to  (subscript  L  for  “lower 
bound”),  which  we  define  as: 

cL(t)  =  e0H(t),  (4.3.12) 

where  e0  is  the  unique  value  to  which  the  solution  jumps  at  t  =  0  (if  it  does 
so),  given  in  equation  (3.3.4),  and  f{  {  t)  is  the  Heaviside  step  function.  The  next 
lemma  validates  Ej  as  a  lower  bound. 
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Lemma  4.3.7.  If  £in  vanishes  for  t  <  0,  is  monotone,  and  satisfies  £jn(r)  >  £p(r) 
for  all  t  <  t,  then  e  0  u  t  (  t)  >  e  p(  t),  for  t  >  0,  when  the  material  is  not  purely  elastic. 


Proof:  From  Lemma  4.3.2  we  have 


{J'*e  i„)(<)  >  {J1  *eL)(t), 


(4.3.13) 


with  equality  if  and  only  if  sin  and  are  identical.  Then,  since  convolution  with 
the  Heaviside  step  function  is  just  an  integration: 


(J1  *eiD)(t)  >  J(t)ec 


>  J  0£  Q  11  (t). 


(4.3.14) 


for  materials  which  are  not  purely  elastic.  Thus, 

fout(')  =  f(u-(J  1  *  ein)(t))  >  F(u-J0e0n  (0)  =  £t(0,  (4.3.15) 

from  the  definitions  of  £out,  £<,>  aad  eL<  an<^  the  m onotonicity  of  F.  C 

Lastly,  for  the  output  to  be  greater  than  the  input,  we  have 

Lemma  4.3.8.  If  ein  vanishes  for  t  <  0,  is  monotone,  and  not  identically  zero 
for  I  >  0,  and  satisfies  er  j  n  ( r )  <  c  i(r)  for  all  r  <  t,  then  eOut(0  >  £in(0>  " here  eiD 
is  non-zero,  when  the  material  is  not  purely  elastic. 


Proof:  We  use  the  monotonicity  of  /  and  equation  (2.3.8)  to  write: 


/  (ein  )  <  f(£o)  =  ™(£o)£  in-  T  < 


4.3.16) 


with  equality  if  and  only  if  ein  is  identical  to  eL.  Now,  m(e0)  =  U~J0.  and  by 
equation  (2.6.9)  we  have  J0£0//  (  t)  <  (J  1  *  eln)  (  <),  when  the  material  is  not  purely 


elastic.  Therefore, 


/(*•,.( 0)  <  **«.)(«)■ 


(4.3.17 


»n  n  vTyrymT  vr.  :  tv.  vr\prf  \rr’jv  '/yj^vg  ■v.--.-%--.  ’vw.vwj-j-.rv.  m 


Using  the  monotonicity  of  F  and  the  definition  of  tout  we  obtain: 

£in(  0  <  F  (v2{J  '  *  ei»)(  0)  =  £out(  t),  (4.3.18) 


which  is  the  stated  result. 


□ 


4.4  U  >  U0. 

We  construct  sequences  of  functions  using  the  interation  scheme  (4.3.3)  for 
wavespeeds  V  >  U0.  VV e  define  an  upper  bound  sequence  to  be  a  monotone  sequence 
of  functions,  each  member  of  which  is  an  upper  bound  on  all  succeeding  members 
Similarly,  a  lower  bound  sequence  is  one  for  which  each  member  is  a  lower  bound 
on  all  succeeding  terms. 

We  first  construct  an  upper  bound  sequence.  We  fix  a  value  of  /  >  0  and  choose 
an  input  function  sq  which  satisfies  the  conditions  of  Lemma  4.3.5  and  is  continuous 
for  t  >  0.  The  first  iterate  e2  satisfies  e2(<)  <  «q(f)  from  this  lemma.  Further¬ 
more,  this  result  applies  for  all  earlier  times  greater  than  zero;  we  have  e2(r)  <  £q(r) 
for  0  <  t  <  t.  We  recall  that  Eg(t)  >  £/,(<)  =  e0II[t)  for  t  >  0,  with  equal¬ 
ity  only  at  t  =  0.  Therefore,  our  starting  function  £j  satisfies  the  conditions  of 
Lemma  4.3.7,  with  the  result  that  £{,(r)  <  £2(7')  for  0  <  r  <  t.  Moreover,  by 
Lemma  4.3.3,  e2  is  monotone  and  continuous  for  t  >  0,  since  £j  is.  By  continu¬ 
ing  the  iteration  process,  Lemma  4.3.4  guarantees  that  the  succeeding  outputs  are 
less  than  their  associated  inputs.  We  thus  have  a  monotone  sequence  of  monotone 
functions  which  are  continuous  for  t  >  0  and  bounded  below  by  £ t  : 

£/,(r)  <  £„+i(t)  <  £n(r)  <  ...  <  e2(r)  <  «T  (0,  0  <  r  <  t,  (4.4.1) 


with  equality  only  at  r  =  0  and  only  if  £i(0)  =  e„  ( i . e . ,  when  the  starting  function 
takes  the  value  of  the  solution  at  t  =  0).  We  similarly  construct  a  lower  bound 


sequence  by  starting  with  an  input  function  which  is  continuous  for  t  >  0  and 
which  satisfies  the  conditions  of  Lemma  4.3.8,  producing  an  output  e2  w*>ich  is 
larger  than  sq  and,  by  Lemma  4.3.6,  which  is  bounded  above  by  eq  .  The  result 
is  a  monotone  sequence  of  monotone  functions  which  are  continuous  for  t  >  0  and 
bounded  above  by  Eq '■ 

e  X  (r )  <  e2(t)  <  ...  <  en(r)  <  en+,(r)  <  eQ(r),  0  <  t  <  l,  (4.4.2) 

with  equality  only  at  t  =  0  and  only  ifsj(O)  =  £0- 

From  the  discussion  of  section  4.2,  we  have  that  such  upper  or  lower  bound 
sequences  converge  to  limit  functions  which  may,  in  general,  depend  upon  the  start¬ 
ing  functions.  The  following  theorem  guarantees  that  all  such  sequences  generated 
from  our  iteration  scheme  necessarily  converge  to  the  same  function  for  U  >  U0. 

Theorem  4.4.1.  Existence  and  Uniqueness  for  U  >  U0.  If  £i(f)  is  con¬ 
tinuous  for  t  >  0  and  satisfies  the  conditions  on  £in  in  cither  Lemma  4.3.5  or 
Lemma  1.3.8,  then  the  upper  bound  or  lower  bound  sequence  of  functions  produced 
by  the  iteration  scheme  (4.3.3)  converges  to  a  monotone  function  which  is  the  unique 
solution  to  (4.3.1)  among  the  class  of  functions  which  are  continuous  except  for  a 
Unite  jump  discontinuity  at  l  =  0.  Furthermore,  these  sequences  converge  uniformly 
on  the  interval  [0,Z]  for  any  finite  t  >  0. 

Proof:  We  first  remark  that  there  exist  suitable  starting  functions,  Sj,  examples  of 
which  are  e j  and  eq  for  lower  and  upper  bound  sequences,  respectively.  We  prove 
that  an  upper  bound  sequence  converges  to  a  solution.  The  same  arguments  apply 
for  lower  bound  sequences.  We  denote  by  £*  a  function  to  which  any  upper  bound 


sequence  (-4.4.1)  converges.  From  (4.3.3)  we  have: 

e*(0  =  lim  £„+!  (<)=  lim  F  (u  ‘  (  J  '  *  e„)  (  <)) 

=  f(u2  lim  (J  '  *  £„)(/)), 

'  n  — ♦  oo  ' 

since  F  is  continuous.  Then  by  dominated  convergence  [4.2], 

e‘(t)  =  f{u2(J'*  £*)(<)), 

since  the  integrands  are  dominated  by  the  first  one  in  the  sequence,  and  this  first  one 
was  chosen  to  be  integrable.  (For  lower  bound  sequences,  Lemma  4.3.6  provides  the 
bound.)  Thus  c*  is  a  solution.  This  limit  function  is  monotone  since  it  is  the  limit 
of  a  convergent  sequence  of  monotone  functions  [4.2].  Now,  a  monotone  function 
can  be  discontinuous  only  if  it  has  finite  jump  discontinuities  [4.2].  Since  e*  is 
a  solution,  it  vanishes  for  t  <  0  and  can  have  only  the  single  jump  discontinuity 
which  occurs  at  t  =  0  (Property  5,  §§3.3).  Therefore,  the  monotonicity  of  £  *  implies 
its  continuity  for  l  >  0. 

We  now  prove  that  the  solution  to  (4.3.1)  is  unique  among  functions  which 
are  continuous  except  for  a  finite  jump  discontinuity  at  t  =  0.  We  will  show  that 
an  assumption  of  two  different  solutions  leads  to  a  contradiction.  We  begin  by 
considering  the  difference  [e: ^  —  £o\(  t)  of  any  two  solutions,  identical  or  not,  at  a 
time  t  which  is  larger  than  the  greatest  time,  say  t0,  before  which  the  solutions 
agree  (all  solutions  at  these  wavespeeds  agree  until  at  least  t  =  0+): 

1*1  -  co|(  I)  =  |  f(u2{j'  *£,)(<))  -  f(u2{J  '*  £2)(<))  j- 

Using  equation  (2.3.16)  to  bound  the  right  side,  we  obtain: 

\£y  -  e2\(t)  <  M  (u2(J  '  *  £p)(t))u2\(J  '  *  (£}  -  £2))(t)r,  p  e  {  1,2}, 

with  equality  if  and  only  if  the  convolutions  are  equal  at  t.  Such  an  equality  requires 
either  that  £\(r)  is  identical  to  £2(r)  f°r  T  <  F  or  that  the  difference  (e]  -  f2)(r)  *s 


not  of  one  sign,  where  it  is  non-zero,  for  r  <  t,  in  view  of  Lemma  4.3.2.  With  the 


use  of  a  triangle  inequality,  we  pass  the  absolute  value  inside  the  integral  to  get: 

k i  -  e2l  10  <  A/(t/2(J  '  *  ep)(t))u2(J  '  *  \et  -  e2\ )(/), 

with  equality,  now,  if  and  only  if  £i(t)  is  identical  to  £2(t)  for  t  <  t . 

We  now  assume  that  sq  differs  from  e2  somewhere  on  [<<>,<].  The  continuity 
of  e  |  and  e2  ensures  that  the  difference  |cj  —  e2\  achieves  its  maximum  somewhere 
on  [/„,/].  By  hypothesis,  e  l  (t0)  =  e2[t0)\  hence,  j  —  e2\  is  not  equal  to  its  maximum 
on  the  entire  interval,  and  we  have  the  strict  inequality: 

lel_e2l(0  <  M  (  U2  {J  '  *  £  a)  (  t])  U2  J  {t  -  1 0)  ill  ax  |e  i  —  |  ( r ) .  (4.4.3) 

V  '  rel(0,t) 

In  obtaining  (4.4.3),  we  have  also  used  the  monotonicity  of  J  and  the  fact  that  £j  —  e2 
vanishes  for  t  <  t0.  Now,  is  a  solution  by  hypothesis.  We  have  from  equa¬ 
tions  (3.1.10)  and  (2.3.13): 

M  (U 2  J  1  *  e  p)  =  M(f(ep))  =  -L-  <  (4.4.4) 

7ri  l  ^  $  )  U  Jo 

The  inequality  in  (4.4.4)  is  obtained  for  increasing  m  since  ep(t)  >  £0  f°r  *  >  0 

o 

and  rn(e0)  =  U  "  J 0.  Thus, 

M  (u2(J  '  *  ep)(  <))  U2Jo  <  1  for  all  t  >  0. 

Since  ./ (  /)  is  continuous  for  t  >  0,  there  is  a  time  t  j  >  t0  such  that 

M  (U2{J  '  *  efi)(tl))u2J  («!-/,)  <  1.  (4.4  5) 

Using  this  in  equation  (4.4.3),  we  have 

kl  “  ^2 1( f  1 )  <  max  lel  ~  (4  4.6) 

re  (<o,r ,  I 
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Hence,  the  maximum  occurs  at  a  time  <j  <  tj.  Equation  (4. -4.5)  is  also  satisfied 


at  t  j  and  we  have: 


kl  -  £2 1  ( ^  i )  <  \£\  ~  e  2 1  ( ^  i )  • 


This  contradiction  follows  directly  from  the  earlier  assumption  that  £j  /  e2  some 
where  on  \t0,t}.  We  conclude  that  the  solution  is  unique  in  the  class  of  functions 
considered. 

Lastly,  the  uniform  convergence  on  an  interval  [0,  <)  of  the  sequence  of  iteration 
functions  follows  from  the  continuity  of  e*  and  each  member  of  {en}°°_,,  according 
to  the  following  version  [4.2]  of 


Dini’s  Theorem.  Let  g,  J\,  be  continuous  functions  on  [a,  6] 

such  that  f  i  >  f  2  >  ...  and  /„  ( t )  — *  ?(f)  as  n  — *  oo  for  all  i  G  (a ,  6). 
Then  /„  — ►  g  uniformly. 


W'e  note  that  Dini’s  theorem  also  applies  for  lower  bound  sequences  {/n}^°=1  by 
considering  the  upper  bound  sequence  whose  members  are  gn  =  g  —  /„  converg¬ 
ing  to  zero  in  the  statement  of  the  theorem  above.  This  concludes  the  proof  of 
Theorem  4.4.1.  C 


4.5  Ue  <  V  <  U0. 


Recall  in  Section  3.3,  we  chose  as  a  candidate  solution  for  these  values  of 
wavespeed  a  function  which,  for  t  — *  oo  ,  is  asymptotic  to  an  exponential  solution 
of  the  linear  problem  in  which  /(c)  =  c.  With  this  in  mind  and,  since  solutions  for 
these  values  of  V  have  no  finite  jump  discontinuities  (§§3.3),  we  consider  continuous 


solutions  of  the  form: 


e(  t)  =  w  (  f )e  , 


(4  5  1 


where  tv  is  continuous  and  satisfies 


I  ini  w{t)=  1  ,  (4  5  2) 

(-»  -  oo 

and  r  is  the  unique  solution  to  (3.2.3).  For  simplicity,  and  without  loss  of  general¬ 
ity,  we  have  set  the  time  shift  tc  in  equation  (3.2.2)  equal  to  zero.  We  again  prove 
that  upper  and  lower  bound  sequences  from  the  iteration  scheme  (4.3.3)  converge  to 
a  solution  of  (4.3.1).  We  will  also  prove  that  the  solution  is  unique  among  functions 
which  have  the  form  of  (4.5.1)  for  these  wavospeeds. 

We  begin  with  the  proofs  of  some  preparatory  lemmas,  the  fi  rst  of  w  h  ich  shows 
that  crt  (which  solves  the  linear  problem)  does  not  solve  the  nonlinear  problem. 
Rather,  it  starts  a  decreasing  (upper  bound)  sequence  when  fj  =  ert  in  the  iteration 
scheme.  Additionally,  it  is  an  upper  bound  for  increasing  sequences  when  <  ert . 

Lemma  4.5.1.  If  £;n(r)  <  eTT  for  all  r  <  t,  where  r  satisfies  (3.2.3),  thee 
£out(  0  <  erl 

Proof:  From  Lemma  4.3.2  we  have 


with  equality  if  and  only  if  ein 

O 

u-j 


J'  *  ein  <  J  '  *  e 

is  identical  to  erl . 

I  ft  rr  2  I  /  V 

*  e  =  U  r  J  (  r )  e 


rt 

Now 

rt 

—  e 


rt 


sinro  V"rJ(r)  -  1  Therefore, 

£« .n(')  =  r(u-(j  1  *  e in ) (  0 )  <  r(cr<)  <  er‘. 

wh'Ti'  \v  <•  have  used  /•’(/)  <  /  to  obtain  the  last  inequality  C 

To  obtain  a  lower  bound  for  upper  bound  sequences  (and  at  the  same  time 


a  function  which  starts  a  lower  bound  sequence),  we  use  the  lower  bound  on  the 


function  F  given  in  equation  (2.3.19)  and  repeated  here  for  convenience: 


F(f)  >  f  -  kf 1 ,  k  >  0,  7  >  1  (-1.5.3) 

fly  using  this  bound  when  ein  =  ert  in  (4.3.2),  we  find  that 

£out(0  >  «ri  -  kU"irJ  (qr)e7r' . 

VV e  note  that 

U2  ^rJ  (7r)  <  1  for  7  >  1,  (4.5.4) 

when  U'rJ(r)  —  1  because  rj(r)  is  completely  monotone  (cf.  equation  (2.5.5) 

and  §§3.2).  To  obtain  a  lower  bound  for  more  general  inputs,  we  consider  input 
functions  at  least  as  large  as  eb  which  we  define  as: 

MO  =  “  Ce1Tl,  7  >  1,  (4.5.5) 

for  some  positive  constant  C .  If  we  choose  C  to  be  the  positive  constant  given  by 

k 


C  = 


(4.5.6) 


1  -  l/27rJ(7r)  ’ 

then  is  the  desired  lower  bound,  at  least  on  some  interval  (— oo  ,<),  according  to 


Lemma  4.5.2.  If  ein(r)  >  et , ( r )  for  a II  t  <  t,  where  r  satisfies  (3.2.3)  and  C  is 
given  in  (4.5.6),  then  there  exists  a  finite  tb  such  that  eDut  (  0  >  eb  (  0  f°r  a/i  1  <  t 1,  ■ 


Proof: 

«out(0  >  f(u2(J'  *  eb)(t)) 

—  F  (^ert  -  CU  "  7r  J  (7r)  c7r(  ^ 

>  crt  -  C  U  “  7r  J  (7r)eirl  —  k‘i’"fr([l  -  CU  ~~jrJ  ( 7  r ) « ^  "*  ’*rl]7  , 

where  we  have  used  (4.5.3)  to  bound  F  and  have  factored  clrl  out  of  the  bracketed 

quantity.  'Phis  remaining  bracketed  quantity  raised  to  any  power  7  >  1  is  a  well- 

defined  positive  number  smaller  than  unity  for  t  <  tb.  w  here 

In  (C  U  -  irJ  (  qr) ) 

tb  - - .  (4.5.7) 

(7-1 )r 
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We  note  that  tm  <  tb,  since  equation  (4.5.4)  and  7  >  1  imply  that  the  argument  of 
the  logarithm  in  (4.5.9)  is  larger  than  that  in  (4.5.7).  We  now  assert  that  e L  is  the 
desired  lower  bound  for  upper  bound  sequences  and  that  it  starts  a  lower  bound 
sequence  upon  iteration. 


Lemma  4.5.3.  If  e  in  (  t  )  >  e^fr)  for  all  t  <  t  and  e 1  defined  in  (4.5.8),  then 

£<>ut(  t)  >  £[A  0- 


Proof:  If  ein  =  then  £out(f)  >  £;_(/)  for  all  I  <  tm  by  Lemma  4.5.2.  Fur¬ 

thermore,  cout  is  monotone  and  continuous  according  to  Lemma  4.3.3,  since  s  ^  is 
both  monotone  and  continuous  by  construction.  Therefore,  £out(()  >  e  t  (  0  f°r 

t  >  >,n  If  £  jn  is  not  identical  to  cj  ,  the  result  is  immediate  from  what  we  have  just 

proven  and  Lemma  4.3.4.  □ 

'I'he  foregoing  lemmas  of  this  section  provide  starting  functions  and  bounds 
for  iteration,  all  of  which  have  the  desired  properties  (4.5.1)  and  (4.5.2).  We  now 

have  all  we  need  to  prove  the  main  theorem  for  these  wavespeeds. 


Theorem  4.5.1.  Existence  and  Uniqueness  for  U e  <  U  <  U0.  If,  for  all  real  I. 


£  1  [  t)  —  crt  w  here  r  satisfies  (3.2.3),  or  if  s  j  (  f )  =  £ g(  t)  given  in  (-1.5.8),  then  the 
iteration  scheme  (4.3.3)  produces  an  upper  bound  or  a  lower  bound  sequence,  re¬ 
spectively,  which  converges  to  a  continuous  monotone  function  on  (—  oo,  +  oo).  In 
either  case,  the  limit  function  is  the  unique  solution  to  (4.3.1)  among  the  class  of 
continuous  functions  which  are  asymptotic  to  ert  for  t  — *  —  oc  .  Moreover,  these 
sequences  converge  uniformly  on  the  interval  (— oo,f]  for  any  finite  t. 


Proof:  If  £ i  =  ert  (ej  =  s the  use  of  Lemma  4.5.1  (Lemma  4.5.3)  followed  by 
repeated  application  of  Lemma  4.3.4  shows  that  the  iteration  scheme  produces  a 
decreasing  (increasing)  sequence  of  functions  bounded  below  by  e ^  (above  by  er'  \ 
according  to  Lemma  4.5.3  (Lemma  4.5.1).  Each  member  of  the  sequence  is  mono¬ 
tone  and  continuous  according  to  Lemma  4.3.3,  since  ej  is  monotone  and  continu¬ 
ous.  Thus,  the  sequence  converges  to  a  monotone  function  e*  on  (  —  oo,+oc).  This 
limit  function  is  a  solution  to  (4.3.1)  by  dominated  convergence,  as  in  the  proof 
of  Theorem  4.4.1.  Since  a  solution  at  these  wavespeeds  cannot  have  any  finite 
jump  discontinuities  and  since  a  monotone  function  can  be  discontinuous  only  if 
it  has  such  discontinuities,  we  conclude  that  e*  is  continuous.  Furthermore,  e’  is 
asymptotic  to  ert  as  t  — >  —  oo  since 


er‘  _  Ce7rl  <£*(/)<  ert,  t<l 


where  ('  is  given  in  (1.5.6)  and  lb  in  (4.5.7).  Therefore, 


1 1  m - —  1 . 

t  — *  —  oo  er 


To  prove  that  e*  is  the  unique  continuous  function  having  this  asymptotic 
property,  we  write  e*  in  the  form  given  in  (4  5.1)  and  (4.5.2)  to  obtain  a  bound  on 


any  such  convolution  U~J  *  e  : 

/oo 

<r*(t  -  7)  dJ  ( r ) 

0 

f  °° 

=  er  1 U  “  v>  *  ( I  —  t  )  e  rr  dJ  (  r  ) 

J  n 

<  ert  in  ax  |  w  *  (  t  )  \  V  “ 

7  <  < 

=  e  r  1  ni  a  x  |  w  *  (  t  )  | ,  (4.5.10) 

r  <  t 

since  the  last  integral  is  rj(r )  and  U~rJ(r )  =  I  by  the  choice  of  r.  Now.  we 
assume  that  there  are  two  different  solutions  tj  and  s  and  obtain  a  contradiction. 
We  consider  the  difference  1^  —  e2|(f)  as  in  the  proof  of  Theorem  4.4.1;  using  the 
bound  in  (4.5. 10)  w  e  obtain: 

Im  -  £j|(  0  <  crt  M  ( U2(J  '  *  c  ^  )  (  <))  niax  |«',  -  w2|(r),  0  €  {1,2}, 

with  strict  inequality  since  uij  and  u’2  agree  at  -  oo  .  When  we  divide  by  ert  and 
use  the  fact  that  M  decreases  from  A/  (0)  =  1,  we  find  that  tuj  ~  tr2  must  satisfy 

|tc .  -  tii  o  |  (  / )  <  m  ax  |  tit  j  -  «io|(7)  for  all  t. 

T<  l 

From  arguments  identical  to  those  following  equation  (4.4.6)  in  the  proof  of  Theo¬ 
rem  4.4.1,  we  find  that  the  assumption  of  two  different  continuous  solutions  leads 
directly  to  a  contradiction. 

To  prove  that  the  sequence  { c  n  }  ^_  j  converges  uniformly  to  e*  on  any  in¬ 
terval  (  oc,/|.  we  note  that  |e„  -  t*|(r)  <  trr  since  ( r)  and  e‘(r)  are  both 

bounded  above  by  rrr  and  below  by  zero  For  any  A  >  (I  and  for  all  ri  >  1  we 

have  (-„  c  *  I  (  7  )  <  A  whenever  r  <  .  where  ~  (In  <s)/r  Thus,  the  convergence 

is  uniform  on  (  x  .  /  A  ]  If  I  >  I  „  the  uniform  convergence  on  |f,v,t|  afforded  by 


j  e~  rT  dJ(r) 


Dun's  Theorem  ((j§4  4)  extends  the  result  to  (-  x-  .  f  | 


r. 


4.6  V  =  U0. 

For  tins  particular  value  of  wavespeed,  our  candidate  solutions  are  continuous 
and  vanish  for  t  <  0.  As  before,  our  plan  is  to  construct  upper  and  lower  bound 
sequences  of  such  functions  and  to  prove  that  they  converge  to  a  solution  of  equa¬ 
tion  (4.3.1).  We  will  prove  that  this  problem  has  a  unique  solution  among  the  class 
of  continuous  functions  which  vanish  for  t  <  0. 

We  begin  by  deriving  an  equivalent  form  of  (4.3.1)  useful  at  this  value  of 
wavespeed.  We  add  and  subtract  J0II  (/ —  r)  inside  the  convolution  in tegral  in  equa¬ 
tion  (3.1.10)  and  evaluate  the  added  term  to  obtain: 

/  (^)  =  U  J  0e  +  U  (J  —  J  o)  *  £  •  (4.6.1) 

Using  V  ~  J  „  =  1.  we  write  this  in  terms  of  the  reduced  strain  curve  {(e)  -  e: 

/(£)  -  e  =  X  '  *  f .  (4  6  2) 

w  h  i>  re 

J  (t)  -  J o 

\  (  < )  =  - - - //(f)  (4.6.3) 

J  0 

is  the  normalized  compliance  which  vanishes  for  l  <  0;  \  is  identically  zero  if  and 
only  if  the  material  is  purely  elastic.  We  note  that  the  reduced  strain  curve  is 
monotone  and.  therefore,  lias  an  inverse,  since  {(e)  —  e  =  (rn(e)  -  l)e,  each  factor 
of  which  is  strictly  increasing  in  e.  We  assume,  as  in  equation  (2.3.18),  that  there 
are  constants  k  >  0  and  >  I  such  that 

f(e)  -  e  -  ke’’ ,  e  (  0  (4  6  4) 

In  terms  of  the  scaled  strain  u(  1)  defined  by 

i 

u(  f)  =  k  ■>-'  e  (  1 1,  (465) 


cws 


>  ,  ' 


we  consider  the  problem  equivalent  to  (-4.3.1)  given  by 

«(0=G((X',*«)(0).  (-4.6.6) 

where  C!  is  not  the  modulus;  rather,  we  now  denote  by  G  the  strictly  increasine 
inverse  of  the  reduced  strain  curve  after  scaling: 

G(g)  =  glh  L(g),  1(0)  =1,  (-1.6.7) 

in  which  /,  is  an  otherwise  unspecified  slowly-varying  function.  Solutions  u(f) 
to  (-1.6.6)  with  (-1.6.4)  and  (4.6.5)  are  solutions  to  (4.3.1).  The  corresponding  iter¬ 
ation  scheme  is: 

«h+i(  0  =  G((X  '  *  “«)(  <))  i  «  =  1,2,  ...  .  (4.6.5) 

1 _ 

The  quasi-elastic  solution  Uq  =  k^~leQ  satisfies 

«<?(<)  =  C(x(0«<?(0)-  (4.6.9) 

We  continue  with  our  now  familiar  attack  on  the  problem.  The  iteration 
scheme  (4.6.8)  is  amenable  to  a  set  of  statements  analogous  to  Lemmas  4.3.3,  4.3.4. 
and  4.3.6  which  allows  us  to  assert,  without  a  redundant  proof,  that  Uq  starts 
a  decreasing  sequence  of  functions  upon  iteration  and  that  it  may  be  used  as  an 
upper  bound  for  an  increasing  sequence. 

Lem  m  a  4.6.1.  If  u  (  is  monotone  and  continuous,  vanishes  for  t  <  0,  and  satis¬ 
fies  11,(7)  <  Uq(t)  for  all  t  <  I.  then  ue(l)  is  a  continuous  monotone  function  of  t 
and,  when  the  material  is  not  purely  elastic,  u  ■_>  satisfies  „  _>  (  l )  <  ti  q  (  /).  for  /  >  0 

We  need  only  exhibit  a  function  U/  which  is  a  lower  bound  for  a  sequence 
started  by  u  q  and  which  starts  an  increasing  sequence  bounded  above  by  uq.  To 


do  so,  we  use  a  lower  bound  on  the  normalized  compliance  Since  x(0)  =  0  and 


its  derivative  is  completely  monotone  (§§2.5),  x  is  convex  down.  For  small  enough 
times,  it  is  bounded  below  by  a  line  of  slope  a: 

X(  t)  >  at  for  0  <  t  <  ta,  (4.6.10) 


w  here 


0  <  a  =  ^—2^-  <  x'(0+)  <  oo  . 


t  a 


(4.6. 11 ) 


We  also  use  a  lower  bound  on  the  function  G  given  in  (4.6.7).  Since  L(g )  — •  1 

as  g  |  0,  there  are  positive  constants  6  and  gs  such  that 


G(g)  >  £? 1  ^  ( 1  -  5)  for  all  g  <  gs. 


(4.6.12) 


Furthermore,  we  can  always  take  6  <  1  since  G  increases  from  C(0)  =  0.  Using  these 
bounds  on  \  and  G ,  we  will  show  that  a  lower  bound  for  decreasing  sequences  (and 
a  starting  function  for  an  increasing  sequence)  is  given  by 

r  _l_ 


( o  =  ■( 


rf~', 


t  <  tr\ 


(46.13) 


l  rt?  ‘  ,  t  >  tr; 


for  small  enough  values  of  the  positive  constants  r  and  tr. 


Lemma  4.6,2.  There  are  positive  constants  r  and  tr  such  that,  if  tij  is  a  continuous 
monotone  function  which  van  is  hes  for  t  <  0  and  w  hicb  satisfies  u  i  (r )  >  u  ^  ( r )  for  all 
r  <  t.  with  u  j  defined  in  equation  (4.6.13),  then  uo(0  ,s  a  continuous  monotone 
function  oft  when  the  m  aterial  is  not  purely  elastic,  and  «2  (  t)  >  u  i(  t)  for  all  t  >  0 
Til  rt  her  in  ore,  u  q  (  t )  >  u  /  (  t)  for  all  I  >  0 . 


Proof:  Fix  A]  6  (0,1)  and  lot 


(1  -  «,)’-• 
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Wo  choose  a  value  of  r  <  and  define 


u  r  (  < )  =  r  t 1~l . 

Using  the  bound  on  \  given  in  (4.6.10),  we  have 

7  —  1  1 

(x  '  *  «r)  (  0  >  ar  - <■»-',  0  <  t  <  ta. 

7 

From  the  lower  bound  on  G  in  (4.6.12),  there  is  a  time  >  0  such  that 

G  (  ( x  '  *  «  r )  (  t ) )  >  (  a  r  - j  t  1  ( 1  -  6  j ) ,  0  <  <  <  m  in  {  /  j . 

>  «r(  <), 

since  r  <  r  ( . 

Starting  anew,  we  fix  (5q  €  (0,1)  and  define 

rQ  =  (1  -  5q)t-'  . 

Now  let  ur  be  as  in  (4.6.14)  with  some  value  of  r  <  rg  Recall  that  t/^  = 
Using  the  lower  bound  on  \  and  the  monotonicity  of  G,  we  have 

uQ(t)  >  G  (at  uQ(t)),  0  <  t  <  ta. 

The  lower  bound  for  G  ensures  there  is  a  tq  >  0  so  that 

«q(  0  >  (at  Uq(1))1I~i  { 1  -  6q),  0  <  t  <  min  {tq,ta}. 

Solving  for  Uq  , 

uQ(t)  >  (af)^(l  -  6q)^ 

>  «r(  n, 

sin  co  r  <  rg  . 

\\  ith  the  above  results,  we  have  conditions  which  allow  us  to  ensure  that 
and  G(\  '  *  nr)  >  « r .  We  now  pick  values  of  r  and  tr  which  satisfy 


(4.6.14) 


t  a  } 


G  ( \  u  q  )  • 


Uq  >  U  ■ 


and  define  uL(t)  according  to  (4.6.13).  If  «j  is  identical  to  u  ^  we  have  just  proven 
that  u,  =  G  ( x  '  *  «i)  >  «i  for  t  6  (0,  f  r ].  Since  u  i  is  monotone  and  continuous  by 
construction,  no  's>  loo,  and  the  result  is  obtained  for  all  l  >  0.  For  functions  ti  j  not 
identical  to  u  j  ,  the  result  follows  from  what  we  have  just  proven  and  a  statement 
for  this  iteration  scheme  analogous  to  Lemma  4.3.4.  Similarly,  for  these  values  of  r 
and  tT  we  have  uq  >  u  i  for  t  €  (0,<r].  Since  uq  is  continuously  increasing,  the 
result  holds  for  all  t  >  0.  □ 

We  are  now  prepared  to  state  and  prove  the  main  theorem  for  this  wavespeed. 

Theorem  4.6.1.  Existence  and  Uniqueness  for  U  =  U  0.  If  u  i(  t)  satisfies  the 
hypotheses  of  either  Lemma  4.6.1  or  Lemma  4.6.2,  then  the  iteration  scheme  (4.6.8) 
produces  an  upper  bound  or  lower  bound  sequence,  respectively,  which  converges  to 
a  monotone  function  satisfying  equation  (4.6.6)  for  the  scaled  strain  u.  The  corre¬ 
spond  in  g  strain  e  from  equation  (4.6.5)  is  the  unique  solution  to  (4.3.1)  among  the 
class  of  continuous  functions  which  vanish  for  l  <  0.  Furthermore,  these  sequences 
converge  uniformly  on  the  interval  [0,fj  for  any  finite  t  >  0. 

Proof.  The  use  of  Lemma  4.6.1  (Lemma  4.6.2)  followed  by  repeated  application  of 
a  statement  for  the  iteration  scheme  (4.6.8)  analogous  to  Lemma  4.3.4  produces  a 
decreasing  (increasing)  sequence  of  functions,  each  member  of  which  is  monotone 
and  continuous.  The  sequence  is  bounded  below  by  it i  (above  by  Uq )  and,  there¬ 
fore,  converges  to  a  monotone  function  u*  on  (— oo,+oo).  By  dominated  conver¬ 
gence,  it*  is  a  solution  to  (4.6.6);  the  cor  res  ponding  strain  e*  satisfies  equation  (4.3.1) 
and  is.  therefore,  continuous  according  to  the  arguments  in  Theorem  4.4.1. 

Uniqueness  of  the  solution  is  obtained  from  the  original  formulation  (4.3.1) 
in  terms  of  F,  the  proof  is  identical  to  that  given  in  Theorem  4.4.1.  In  this 
case,  we  note  that  equation  (4.4.4)  holds,  since  Cp(t)  >  e0  =  0  for  t  >  0  and 
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since  m(0)  =  1  =  U“J0. 


Lastly,  the  u  n  iform  convergence  of  the  sequence  of  functions  on  an  interval  [0.  t ' 
follows  from  Dini’s  Theorem,  as  in  the  proof  of  Theorem  4.1.1.  [3 

4.7  A  theorem  for  bounds  via  trial  and  error. 

In  the  previous  sections,  we  have  exhibited  means  for  constructing  solutions 
from  sequences  of  functions  which  globally  bound  the  solution,  either  from  above  or 
below.  The  first  members  of  these  sequences  provide  the  loosest  bounds.  Moreover, 
the  lower  bounds  we  used,  e ^  and  u do  not  resemble  the  global  behaviour  of  the 
solutions.  For  the  purpose  of  applications,  bounds  for  the  solution  up  to  a  time  T 
can  be  useful  if  they  can  be  made  tight  to  some  degree  without  requiring  many 
iterations  in  the  constructive  iteration  scheme.  The  following  theorem  asserts  the 
usefulness  of  a  trial  and  error  approach  to  finding  starting  functions  which  provide 
tighter  bounds  on  the  solution  upon  iteration.  It  allows  us  to  make  a  guess  for  the 
starting  function  on  an  interval  (-oo,Tj]  and  test  its  iteration.  If  the  iterate  moves 
up  (down)  on  the  sub-interval  ( —  oo  ,  T ),  for  T  <  T\,  then  the  iterate  is  a  lower  bound 
(an  upper  bound)  on  the  solution  until  T.  In  what  follows,  we  denote  by  x(t)  a 
generic  solution,  either  the  strain  e{  t)  or  the  scaled  strain  u(t)  and  we  let  <F(rn) 
represent  either  F  (U2  J  1  *  en)  or  G  (\  '  *  u  n  )  as  appropriate  for  the  w  avespeed  under 
con  sid  oration . 

Theorem  4.7.1.  Suppose  x ^  (x^)  is  a  known  upper  bound  (lower  bound)  for  a 
solution  to  (-l.fi.  I )  and  for  which  i’ixy)  <  ly  (‘P(xf)  >  xp).  Choose  a  non-negative 
increasing  function  x\(l)  which  satisfies  Xj(f)  <  Xy(t)  (x\(t)  >  x^{t)j 
for  t  6  (-oo,7’]j  and  which  vanishes  for  t  <  0  if  U  >  U0,  or  which  has  the  form 
of  (3.2.3)  if  Ue  <  U  <  U0.  If  there  is  a  time  T  <  Tx  so  that  one  iteration  produces 


<t>(xi(  /  ))  =  £o(f)  —  £  i  (  0  (x2  <  £  ] )  for  all  t  <  T ,  then  the  eh  osen  fu  notion  £  j  staffs 
a  lower  bound  (an  upper  bound J  sequence  which  converges  to  a  solu  tion  on  ( -  oc  .  T  ] 


Proof:  Consider  the  case  for  an  increasing  (lower  bound)  sequence.  The  first  iter¬ 
ate,  £ 2,  is  bounded  above  by  Xy\ 

*l(0  <  £2(<)  =  $(£,(<)) 

<  $(£,;{())  <  *1/(0.  t  <  T , 

since  xj  <  xo  and  £]  <  Xy  for  t  <  T  <  T\ ,  by  hypothesis.  Repeated  application  of 
a  statement  analogous  to  Lemma  -1.3.4  produces  the  increasing  sequence: 

2- 1  (  0  <  *2(0  <  *3(0  <  ••-  <  £(.;  (  0  •  1  <  T, 

at  times  where  the  functions  are  not  identically  zero.  The  same  argument  applies 
with  reversed  inequalities  for  decreasing  (upper  bound)  sequences  bounded  below 
by  x  1 .  Convergence  to  a  solution  with  the  appropriate  properties  is  proven  from 
the  arguments  in  the  proofs  of  Theorems  4.4.1,  4.5.1  and  4.6.1.  □ 


61 


CHAPTER  5:  STEADY  SHOCK 


AND  ACCELERATION  WAVES 


5.1  Problem  Statement  and  Results. 


We  first  consider,  in  section  5.2,  the  explicit  dependence  of  shock  and  accel¬ 
eration  wave  solutions  on  the  parametric  wavespeed  V.  We  show  that  ;(<;(')  is 
continuous  in  U  for  solutions  which  vanish  for  t  <  0.  In  particular,  the  accelerato  n 
wave  is  the  limit  of  shock  waves  as  U  J  U0.  Along  the  way,  we  show  for  two  shock  or 
acceleration  waves  having  different  wavespeeds  that  the  faster  wave  is  everywhere 
larger  than  the  slower  wave,  after  the  initial  jump.  We  conclude  with  a  theorem 
which  is  based  on  this  result;  it  allows  us  to  construct  a  solution  for  a  wavespeed  L  -., 
from  a  known  solution  at  a  different  wavespeed  Ua ,  when  both  wavespeeds  are  at 
least  as  large  as  the  acceleration  wavespeed  U0. 

In  the  next  three  sections,  we  consider  the  behaviour  for  t  1  0  and  t  — *  oc  of 
solutions  e(  t )  for  shock  and  acceleration  waves;  i.e.,  for  values  of  wavespeed  U  >  L'0. 
We  First  obtain  the  small-time  asymptotic  forms  of  e(t)  -  e0 ,  where  e0  is  the  non¬ 
negative  value  the  solution  takes  at  <  =  0+.  When  the  strain  curve  /  satisfies  the 
assumptions  of  section  2.3,  the  result  obtained  in  section  5.3  for  almost-all  shock 
waves  is: 


e(  <)  -  ~  q  X  (  f),  <10, 


(5.1.1) 


w  here 


f'(e0)  -  [U/Vo)2  ’ 

and  \' (  t)  is  the  normalized  compliance  introduced  in  equation  (4.6.3).  Moreover,  th< 


I u as i-elas t ic  shock  solution  lias  the  saint*  asymptotic  form 


Our  small-time  ri-sult  for  an  am*  lerat  io  n  wave,  for  which  -  <1.  relics  up":, 

the  additional  assumption  that  \(t)  is  regu  larly -vary  in  g  for  I  1  <)  with  pow.-r 

l’{\)  -  ji  e  (n,  |]  We  find  that  t(f)  is  also  regu  la  rly  -  vary  in  p  Its  power  i- 

/’  | ;  )  -  7  €  (0,  oo  ) .  w  here 


in  which  is  the  dominant  nonlinear  power  in  the  reduced  strain  nine  /I-:)  -  f.  : 

small  values  of  as  in  equation  (40.4)  1 11  terms  of  the  scaled  sti  ai.’i  11  I  I )  define’, 

in  equation  (I  (>.51.  the  result  we  obtain  in  section  5  I  for  accelerat  1  *  *  ti  waves  is 


«( t)  ~~  (it  \  ( /))  .  I  |  0. 


w  h  e  r  e 


(P  +  Q)' 

Here  the  quasi-elastic  solution  differs.  The  result  is 


«q(  t)  ~  (\(  t))  .  (  |  0 


(5  1  5  | 


(5.1.6) 


The  final  section,  5.5,  is  devoted  to  examples  of  shock  and  acceleration  wave? 
constructed  for  specific  materials  from  a  numerical  approximation  to  the  iteration 

scheme. 


v:-w.*K-:v 


5.2  Dependence  of  solutions  on  the  waves  peed  I  ■ 


\\  c  ii  ii  s  id  <■  r  I  In-  <■  v  |  >  1 1 «'  1 1  dependence  of  solutions  -  ■  ■  n  the  pa  ra  in  e  (  re'  wave- 

s  p . I  /  a  II  it  we  write  c  |  /  .  f  ’  )  We  sometimes  use  the  »<  1 1  a  !  i  -  ,  (  <1  -It  I  ,]  | 

We  have  shown  in  Properties  1  and  5  of  section  t  d  that  t  Ii  •  equilibrium  v  a  I  u  • 

and,  for  shock  and  a  r  ce  |era  t  io  n  waves,  the  initial  jump  value  are  stiictlv  in  c  i  i-asin  e 

functions  of  th*’  wai<'speed  It  is  natural  to  ask  if  such  a  statement  i-  true  f.'t 

s  h  o  e  k  and  a  e  c  e  |e  r  a  t  n  ■  n  waves  f..  i  all  times  where  the  solution  is  n  ■ .  n  •  /  •  r  o  Indeed  i  • 

Is  a  i  ■  o  id  I  II  1'  to 

Lemma  5.2.1.  Among  solutions  w  huh  vanish  nl <•  ii  t  n  .i  1 1\  [nr  I  •  d.  if  I  ,,  <  I  , 

Ihm  -|  l.  r  .)  <  It  t  ,,)  for  a  II  t  >  0 

Proof  Solutions  which  vanish  for  /  <  <)  are  either  shock  or  ac<  <-|e  ra  t  ion  waves  and 
correspond  to  wavespeeds  I  >  For  such  solutions,  the  difference  c„(()  -  t  ( ,(f 

is  continuous  for  t  >  0.  according  to  Theorems  4.4  I  and  1  6  I  Since  ea(0)  <  f(,(0| 
for  „  <  f'j,  (jj§3  .  the  difference  ea{  t)  -  £(,(  t)  is  either  negative  for  all  t  >  0.  and 

the  proof  is  complete,  or  it  is  zero  at  some  finite  time  t0  >  0.  We  will  show  that 
the  latter  case  is  contradictory.  We  subtract  the  equations  satisfied  by  f „  and  f (, 

ea(l)  -  e6(l)  =  f(u‘(J  '  .  ea)[t))  -  F  (  V*  ( J  '  •  e6 )  (  t) ) 

<  F(U‘(J  £„)(<))  -  *  et)(l)). 

since  F  is  strictly  increasing  and  V a  <  f/j,,  by  hypothesis  I  sing  equation  (2  3  Mm 
on  the  rig  lit  side,  we  obtain: 

ea(t)  ~  et(0  <  '*  *«)(/))  V;  (J  '  •  -  **)><  » 

Now  ,  assu in  e  that  t0  is  the  first  tint  e  for  w  h  ich  e„  -  eb  is  zero  Since  r„  (  r  )  -  e 6(  r )  <  ti 
for  allO  <  t  <  t„,  the  convolution  on  the  right  side  is  negative  Furthermore,  M  i> 
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>>>: 


p  •  >Ml  IV  .■  a  ii  .1  wo  li  ;t  v  e  at  I  , 


ii  <  i  M  (  I  ,  M<  I )  1 l  •/  '  •  i  ■  ,  --HIM 

N\  h  11-  )l  ,  ,1 II  1  l  ;t<l  |C  ts  t  ll  <•  .isSU  III  p  t  loll  r  fl  1  t  -  ,  ( t  I  -tit  I  " 

\S  c  Ii  .  i  w  l|  SI'  t  ll  I-  results  of  t  he  |>  re.  -Ill  III  I’  l<  -  III  III  ;i  I  ,  | ,  T  ■  •  V  « ■  I  J !  lit  -  I  I  I  )  I  -  "11 

1 1  ii  ii  •  i  M  ~~  in  I  under  r  <  •  ->  t  r  i  tod  ,  ir<  ii  hi --t  .i  ii  • 

T  Ii  e  c»  r  p  m  5.2.1.  ,A  m  ,>  /i  g  s ,  >  l  u  ( i*  >  11  s  a  h  i  >  l<  \  a  11 1  *  h  i  *  1  •'  n  1 1  a  I  I  \  f"  r  t  *  it  -  l  t  I  a 

•  ,<||  I  III  null  s  function  I  III  I’ll  rt  ll  II  Ia  r  III''  f ' 'h-1  A  i  I'll  ll  |W  a  t  /  I  'It' 

ll  III  I  I  ' if  S  ll  "  •  k  «  ,111'S  .is  I  i  I 

I'r-.-.f  \\  >•  .ire  again  ‘  o  ll  r  e  l  II  eil  Willi  Wav  es| . 1  S  /  -  I  1  I  I  ■  I  ■  "  II '  II.  •.  It  > 

follow  S  f r,,u I  showing  that  -  |  f  /  )  i'  •  > •  > t  h  r ig Ii i  •  o n  i  i n  u ■ . u  s  a  n  ,1  lef i  -c ii t  m  u  o u s  r  ! 

\\  provide  i  he  details  of  the  proof  f,,r  i  ig  h  t  -  c  o  n  t  in  u  it  v  and  in  d  i<  a  t  <  pa  re  n  t  h  >■  1 1  « 1 1  > 

a  in  od  i  lie  a  t  in  n  fo  r  the  similar  p  ro  o  f  ,,f  left  -,  •  o  n  i  in  ii  it  > 

for  a  1 1  \  I  >  I  we  choose  a  convergent  sei  j  ll  e  l|  c  e  of  W  a  4  e  speeds  {I  r  I  ,  _  ! 

for  which  rr,  |  /'  for  each  it .  Theorem  1  t  1  guarantees  the  existence  of  a  'look 

solll  t  loll 

,( I,rn  I  =  /  (  i;ij  1  .  -<  .r„)n  n).  »  1.2.  (■  -  n 

from  Lew  1101  ’>  2  1,  those  monotone  solutions  form  a  decreasing  se<iuonce  of  func¬ 
tions.  ea<li  in  e  m  I,  e  r  of  which  is  hounded  h  e  low  by  -  (  ( .  I  |  1  here  is  a  monotone  func¬ 

tion  7. (f)  such  that  c(  /,/■„)  1  iff)  >  fff.f  )  We  may  pass  to  the  limit  as  I  ..  .  I 
in  equation  (5.2  1),  since  F  is  continuous  and  the  convolution  integrands  are  dom¬ 
inated  by  f(f,f/|)  (for  increasing  sequences  used  111  the  proof  of  |e  ft  -  eon  t  I  ti  II  1 1  > 
as  I '  ,i  ]  f  .  the  integrands  are  dominated  by  the  shock  solution  for  any  wa'espeed 
greater  than  I1  )  In  the  limit,  we  have 

F{  / 1  --  f(/’’2(J  '  .  e)(  n) 

(>S 


lienee  -  (  1 1  is  ,i  shock  so  lu  l  ion  and  li  as  1 1  s  o  n  ly  (i  n  it  e  iliM-unlm  mh  a  I  t  :r  1  i  I.  «■  f.  r  «■ 
a  hi*  It  it  vanishes  id «  n  t  h  a  1 ly  Since  t  is  monotone.  1 1  is  continuous  (, ,  i  /  o  |  i  ,  m 
I  ll  e  II  |l  1  ■  I  I|  e  II  ess  o  f  I  c  ||  so  1 1|  t  loll  s  |)  ro V  id  eel  l)\  'I' ll  eo  fe  III  1  t  I  .  vv  e  ll  a  v  -  I  I  |  -It  I 
i  ii  >1  wo  li  is  o  p  rov  e  n  that 


Inn  f  (  f .  I  ,,  ) 

r„it 


in  -  (  if) 


Ue  similarly  prove  that  there  m  a  fuiietioii  -It)  such  that  ("i  I  •  / 


1 1  m  -■ (  I  I  I  -  -  (  <1  -I  I  I 
f  ■  •  r 


t  I  I  |s  <  ■  .  .  II  I  III  II  o  ||  S  III  / 


I  o  •■l.tain  *  Ii  ••  a  <  <  le  r  a  t  fo  n  wav..  as  the  limit  ..f  shock  w  a  v  .  s  f. ,  i  /  _  f 

(••How  t  li  o  j>  r  ■  i  >f  of  r  ii’  h  t  •  eon  t  in  ii  it  y  a  l.ovo  w  it  h  t  h  e  fi  llo  a  in  it  a  rtfii  mi  on  t  s  ^in.  .  <  a  : 

III  o  II  o  l  . .  ll  .■  III  ■Oil  I  O'  r  of  I  ll  O  SIM  |  U  ell  <  e  {  -  j  (,/„)['  J  vanishes  for  t  i  ( 1  all  <1  s  in  .  I  I, . 

i  i  iii  1 1  s"fui  ion  -  (  1 1  e  a  ii  li  a  v  e  no  d  isro  ii  (  hi  ii  it  ics  at  w  a  cos  peed  / .  wo  have  that  •  I  t )  e 
o  ii  t  i  ii  ii ..  a  s  and  vanishes  for  I  <  U  From  Theorem  I  ii  1  .  such  a  solution  is  unique 


II  o  no  .■ 


in"  <(  i  i •„)  -  ()  -  -i  i  r„) 

'  ,i'  .. 


1  omnia  i  ‘i  I  spawns  another  result  which  is.  perhaps,  more  useful  than  t  h  •• 
l  h  ■  ■  ■  >  r  •  ■  1 1  •  ,t  I  result  of  continuity  of  solutions  in  w  avesp  eed  It  allows  us  to  use  the 
in  foi  in  a  i  ion  i:i  a  known  solution  at  a  particular  wavespeed  to  generate  the  solution 

at  any  other  wavcvp . I.  when  hot  h  w  avspeed  s  in  u  est  ion  are  at  least  as  large  s 

that  "f  ill*-  i)  i  i-  It  r  a  t  io  ii  wave 

Theorem  5.2.2.  Suppose  f  (  f.  >  I’.,  are  given .  with  (  „  /  l  f,  and  /c(  I  t  , 
in  tlii ■  1 1  r  r  a  ( in  ii  srhrnir  ( -t  .'I  .'I  j  If  the  solution  at  wavespeed  I  a  is  know  n  to 
hr  -(  i.l  „  )  -  fair)  for  i  £  (t),f|.  thru  upon  iteration,  f  „  starts  an  i  ipper  hound  fa 


/on  r 

\  rftr 

Pl'-M 

r  *  *  v  •  * 

I  }m* 

n  li  •• 

it  «•  r  . 

I  •  •  1. 

s  in  ( 

i  n  u 


t  r  a  < 

\\  •  * 
in  ( 


r  huun<l)  nf  fumtutns  n  hell  I  ,  *  /  ( /  „  «  /  i  I  I  III s  setjiirnce  (>>n 

#*  s  tin  i  f"  r  m  I  \  <  >  n  |n .  / 1  to  t  he  shuck  <*r  ;t  err  It*  rut  mn  «  a  v  r  so  In  f  /o  //  -  ,  (  ;  )  •  (  ;  .  /  . 

*  f  *■  j»  r <» v i*  t  In*  »■  a s«*  /  * „  >  /  ^  I  li  <•  s a  m  <•  a  ri*  u  iim’  n  I  s  j .  1 1 1  \  with  1  h  »•  1  n  *■  <]  u  a  i  1 '  .  •  ■ 

is.*. |  for  <  f’,,  Suppose  f  „  >  r^,  hii<I  l«*t  •  j  -  in  tin*  1?  ra  l  i«  >  n  sr|i.-i:.» 

n  for  ;  f  [O.f] 

-,(o  mm)  <  r(c,;tJ  '*  -,ho)  -.i.-i  •  1 1  ■■  1 

it-  we  have  used  the  increasing  nature  nf  /'  In  nlit.iiii  the  hi  •' <|  u  a  1 1 1  >  |{  e  |  >  e  i  : .  : 

1 1  h  1 11  produces  a  decreasing  set  u  e  n  ee  w  huh  is  bounded  he  1 1 .  w  |i\  ,|  ;  )  a  •  c . .  t .  i  . 

•  in  ina  t  3  I  and  the  following  ineq  n  a  lit\ 

-  ■(')  r(r?(j '  *  -om)  >  i-  {< '  *  - /. ) < *  i ) 

e  >  -  fi  according  to  Lemma  5.2. 1  when  >  t'h  The  lion  tided  decrea* 
■-e  i|  ii  e  ii  c  e  converges  uniformly  to  the  unique  solution  appropriate  for  the  wave 
•d  f’n  >  according  to  either  Theorem  t  1  I  or  Theorem  l.C>  1 

5.3  Shocks  as  t  J  0 . 

We  begin  |.\  deriving  an  equation  satisfied  by  e(  I )  -  for  /  >  D.  \\  e  sub 
l  /(-.,)  from  each  side  of  equation  (3.1  III] 

/(HM)  /(-  )  c-(-/  '  *  •  ML  d  i 

then  use  (  |  (i  3 )  to  write  ./  (  t )  ./,,[//(()  t  \  (  /  )  ]  in  the  convolution  Integra 

5  3  1)  We  nliUl II 


/ M  0)  -  /(fa)  =  AM  M  +  A(x  '  *  0(0. 


(5.3  2 


.1  II  i  .  I.  (•  '  I  I  ■  I  1  W  .  1 1  . 1  \  i  * 


i  1  '  >  I  t  I 


If  ;  I  '<  \ 


I  '  I 


\  it  •  i  a 


,i/i  if)  II  I  t  I 


II  i  I  >ii  l  li  ■  i  ii  I  •  r  i  i  I  "  f  i  ■  •  ■  I  i  w  .  ■  I  '  n  n 


I  /  I  \  I  f  I  *  (  V  •  ,  (  I  1 


tli,  j  ■  •  .  \  .  i;  r  \s  •  ♦  •  * |  ii 


.tl|"||N  |M  *'<  j  II  it  I  !'  '  II  (  "t  'A  'J  I  V\  r  -‘I't  tin  i  i  ■  ■  i  ?  •  1 1  f1 1  > 


r  i  ,  i  f  *  *  /i  i  Kin  /  ( ■ \  t  n  >  k  \  ,  u 


x  Ii ,  I  ,  w  ||  .1  \  I  I,  .  !  f  I  -  i  A  •  Sinn-  /  i»  assumed  t  •  •  have  a  piecewise  Hillin' 
I .  r  i\  it  iv  ■■  (  |  'a  .  in  av  «  lit'' 


If  '  I  )  A  I  ,  (  I)  I  (  ,  I  n  )  /(-..)  V  I  t)  r  \  (  \  '  •  ,  i|  I 


,s  I, .  i- .  /  is  .i  si.  .u  I  v  -  v  a I  v  ill  p  f  ii  ii  i  t  i'.  ii  fur  w  Ii  i<  Ii  /.(<*)  I  I. '  J  u  at  I*  >  li  I  ’>  •'!  '1  will  I  •  ‘ 

n  ,  I,.  i||i.  result  fii|  .t  1  ill  iwl  .ill  s  |i  i  >  i  K  waves  I  Ii  <•  e  \  r  e  |i  t  1 1  >  II  s  urnir  .1 1  "  a  v  ■■  s  . . I  '  f- 

Ahnli  f  '  |  -  .  |  »  , .  r  f  '  I  ,1  V  cur  res  pun  il  iii  r  .  fn  r  «'\  a  in  |>  I*-  I  "  points  A  ‘  n  '1  I1  ' 

|  I  .■  >|  I  e  |  III  |,  r  a '  l  1 1  •  waves  I  r  a  v  *■  1 1  ii  i’  a  I  ( Ii  ese  e  \  a  ■  t  i.ilucs  "f  «  aViA|H'c  .|  have  :i 

I ,,  ,  I  1 1  .1  1 1  I  1  i  I  V  .if  I.  e  III  l’  e  \  c  1 1  e.l  I II  a  111  a  I  e  III  I  II  "  I  d<>  these  c  a  ses  pi  esc'll  t  a  II  V  l!I  !e|  e-  I  , 

In  a  f  Ii  .  Ill  all'  a  I  |  ■  Il  ■  il  • » 1 1 1  'll  a 

,  now  show  that  the  left  side  of  ei|'iati  n  ( ■<  -1  ~  I  .  'I  n  IV  a  h  ri  t  l>  I  Ii 

left  sole  of  I :i  h|.  IS  ail  HH-reasilin  function  of  ,  and  therefore,  i  n  \  «•  t  t  1 1  ■ 

|ls  derivative  with  respect  l>>  •,  is  |o,sitlve,  a  r  <  c  i  r<i  I  il  g  to  the  relation 


/  |  -  ,  *  1  >  tn  ( -  ,  I  A  .  v\  1 1 1 •  1 1  fc  1 1 1"  <a  f  i .  i  iii  1 1  n  .i  r  i •  •  it  |  'J  A  i,  |  'A  A  ii 

I  A  A  .'I  | ,  aiui  lli''  in  "  n  <>  to  n  H  it  >  cf  in  I  Ii  '■  i  f"  i  '■  «  .•  in  t \  i in  .  r  l  ''ju.iiii.u  i  i  A  'i 

f"  i  ,  (I  I  :i  u  il  will'' 

,  (  n  [<»  \  (  n  i  U  \  1  *  ,  )  (  I  )]  I  (  '■  \  I  M  •  .<  I  \  '  •  ,  I  (  f  I  )  |  i  A  . 

w  hi-ii'  ii  is  Kivi-h  m  t'i|  ii  at  ion  (A  1  2  | . 


/  V  I  \ 

i .  I  /  is  a  s  I"  u  |  v  i;m  i  ii  ('  fu  ii  o  ( i"  n  fo  i  k  /i  e  Ii  /(iii  I 


\  w  ■  "  I!  S  III  ■■  I  t  11  >■  I  llll  It  III  l!  I"'  Il  a  \  I"  II  r  "f  •")  ||  a  I  I"  1 1  (  '.  A  'I  I  f.  !  I  _  II  It'll,  1  .1  ;  : 
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strains.  However,  for  small  q  (large  7),  expansion  of  the  factorials  yields 
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where  (;j!)/  is  the  derivative  of  the  factorial  function  evaluated  at  p  and 
T  =  0.5772  ...  is  Euler’s  constant  [5.1].  The  quasi-elastic  solution  more  closely  ap¬ 
proximates  the  solution  at  small  times  when  the  material  is  weakly  nonlinear  at 


small  strains. 


5.5  Numerical  Examples. 

In  section  5.1  we  presented  asymptotic  results  for  acceleration  waves  which  are 
exact  for  a  material  having  a  specific  strain  curve  and  ;  power-law  compliance.  Such 
a  compliance  grows  without  bound  and  cannot  represent  a  viscoelastic  solid  for  all 
time.  To  better  illustrate  the  global  behaviour  of  acceleration  waves  in  solids,  we 
consider  the  steady  acceleration  wave  problem  for  a  material  having  the  normalized 
exponential  com  pliance 

X(t)  =  Xe  (l  -  Hit),  (5.5.1) 

with  the  strain  curve 

/(e)  =  e  +  ke\  k  >  0,  7  >  1  .  (5.5.2) 

In  terms  of  the  compliance  J  ,  the  equilibrium  value  \  t,  appearing  in  (5.5.1)  is  given 
by  : 

Xe  =  [Je  ~  JoMJo-  (5.5.3) 

We  use  a  simple  numerical  approximation  to  the  iteration  scheme  (4.6.8)  to  generate 
approx  im  a  tion  s  to  a  fi  n  ite  number,  I\  ,  of  terms  be  Ion  ging  to  upper  and  lower  bound 
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sequences  for  the  scaled  strain  of  equation  (-1.6.5).  We  similarly  generate  such  se¬ 
quences  for  shocks  in  m  ate  rials  having  either  a  power-law  or  exponential  compliance 
combined  with  this  strain  curve. 

The  governing  equation  for  the  acceleration  wave  problem  is 

(«(  0)7  =  (x  '  *  «)(  0-  (5.5.4) 

The  iteration  scheme  takes  the  form 

uk  +  1  =  (  ( X  '  *  u  k )  (  1 ) )  ^  i  ^  =  0,1,2,  ...  ,  A'  —  1 .  (5.0.5) 

The  bounding  sequences  are  desired  for  t  G  [0,7’]. 

In  order  to  perform  the  convolution  on  the  iterates  {u  k}  ~  1  ■  we  use  a  piece- 
wise  linear  approximation  to  u  k  and  convolve  this  approximation  exactly  with  \. 
We  partition  the  interval  [0,T]  into  N  subintervals  of  equal  length  given  by 

h  =  T/N.  (5.5.6) 

The  piecewise  linear  approximation  to  a  function  <p  on  (0,7’)  is 

Lj  <p{  t)  =  <Pi  +  -{Pi+i  ~  <P .)  (<-  U)>  ti  <  t  <  ti- n-  (5.5.7) 

n 

w  here 

tf  =  i  h  and  y? ,--£>(< (- )  fo r  j  =  0,1,2,  ...  ,  A' .  (5.5.5) 

The  iteration  scheme  for  the  numerical  approximations  is 

,  =  Lj  (x  '  *  )  ’/7(  0.  A'  =  0,1,2,  ...  .  /v'  —  l ,  (5.5.9) 

«  *,  **  r**  »  ^  (  /)  =  Ijj  u  j.(  t) .  The  upper  bound  and  lower  bound  sequences  we  generate 
in-  tail  i‘d  w  it  h 

un  (  0  =  A  T  u  q  (  t) 


(5.5.10) 


respectively,  where  Uq  is  the  quasi-elastic  solution  given  by 

«<?(<)  =  <e(o,7'l, 


and  u  i  is  t  li  o  lower  bound  from  equation  (4.6.13): 


uL(t)  =  r  t 


t  €  [0,7|. 


(5.5.12) 


(5.5.13) 


For  the  strain  curve  considered,  the  constant  r  is  T]  in  the  proof  of  Lemma  4.6.2: 


(5.5.14) 


with  the  aid  of  equation  (4.6.11). 


From  the  numerical  scheme,  we  generate  the  iterates  ttj;  at  the  mesh  points  t ■. 
Let  !p^  be  any  one  of  the  piecewise  linear  u  ^  .  Its  convolution  at  <,•  is: 


(x  * *  <pN  )  (L)  =  /  <pN 
J  0 


<P  (<i  -  t)  d\(  t) 


=  E  /*. 


where  the  A- - 1 li  integral  in  the  summation  is  given  by 


(5.5.15) 


k  =  /  (<i  -  T)dx(r 

*  A  -  1 

[tk  f  /V  ^  + 

=  J  - 


l(4  ~  r) }  <*x(r), 


(5.5.16) 


in  view  of  (5.5.7)  and  with  the  notation  of  (5.5.8).  We  put  our  results  in  terms 
of  the  non-dimensional  time  t/tT  appearing  in  equation  (5.5.1)  for  the  exponential 
compliance.  With  the  change  of  variables  t/tr  — *  t,  we  have  dx(")  =  Xfe~  dr. 
Using  this  to  integrate  (5.5.16)  exactly,  we  obtain  after  simplification: 
h  =  [('  -  k)'P^k  +  k  P'i-k+  l]  ■  (Xit  -  X  *  -  l ) 

+  ^j-(lP^k+x  ~  '  [(*  +  kk)e~kh  ~  (1  +  (k  ~  l)*)e_(fc~  1)A]  (5.5.1  7) 
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The  numerical  scheme  is  now  completely  defined.  It  consists  of  equations  (5.5.9). 
(5.5.15),  (5.5.17),  and  either  (5.5.10)  or  (5.5.11)  to  get  started. 


We  remark  that  we  chose  this  approximation  to  the  exact  iteration  scheme 
both  for  simplicity  and  since  any  continuous  function  can  br  approximated  arbi¬ 
trarily  closely  by  piecewise  linear  functions  [5.2],  For  a  finite  number  of  terms  gen¬ 
erated  by  this  numerical  scheme,  we  can  gel  arbitrarily  close  to  the  exact  iterates 
by  taking  N  large  enough. 

We  present  the  results  of  this  scheme  as  a  series  of  graphs  of  approximations 
to  acceleration  waves  for  three  values  of  the  nonlinearity  7.  In  Figures  5 . 1  —  5 .  -4 
we  have  plotted  the  normalized  iterates  uk  '  vs.  f/<r-  Figure  5.1  illustrates 

truncated  upper  bound  and  lower  bound  sequences  for  7  =  5/4.  It  shows  iterates 
1-25  of  the  upper  bound  (decreasing)  sequence  started  by  Uq  and  iterates  1-30 
of  the  lower  bound  (increasing)  sequence  generated  by  u  L  .  For  e*ch  curve  in  this 
figure,  201  points  were  used  in  the  computations,  whereas,  101  are  plotted.  The 
large  space  between  the  truncated  sequences  is  the  numerical  approximation  to  a 
bound  within  which  the  acceleration  wave  is  to  be  found.  Our  approximation  to 
this  acceleration  wave  is  illustrated  in  Figure  5.2  along  with  the  quasi-elastic  so¬ 
lution  for  7  =  5/4.  The  approximate  acceleration  wave  is  the  superposition  of 
the  indistinguishable  47-th  members  of  the  upper  and  lower  bound  sequences  of 
the  previous  figure.  Figures  5.3  and  5.4  illustrate  the  quasi-elastic  solutions  and 
similar  approximations  to  the  acceleration  waves  for  7  =  2  and  7  =  5,  respec¬ 
tively.  In  Figures  5.2  5.4,  all  201  points  used  in  the  computation  of  each  curve  are 
plotted.  We  note  that  the  exponential  compliance  is  regu  larly-varv  in  g  as  I  )  0  with 
power  p  —  1 .  The  graphs  for  t  h  ese  t  h  ree  values  of  7  illustrate  app  rox  im  a  lions  for  ac¬ 
celeration  waves  which  are  regularly- varying  for  t  (  0  with  power  q  =  1/(7  -  1);  these 
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Figure  5.2.  The  quasi-elastic  solution  (top  curve)  and  the  superimposed  47-th 
iterates  in  the  upper  and  lower  bound  sequences  representing  the  approximate  nu¬ 
merical  solution  for  the  strain  acceleration  wave  with  7  =  5/4  and  an  exponential 
compliance.  The  ordinate  variable  is  (k / \ e)^  ^ 1  ~  ^e.  The  abscissa  represents  t/tr. 


Figure  5.4.  The  quasi-elastic  solution  (top  curve)  and  the  superimposed  6-th 


iterates  in  the  upper  and  lower  bound  sequences  representing  the  approximate  nu¬ 
merical  solution  for  the  strain  acceleration  wave  with  7  —  5  and  an  exponentia 
compliance.  The  ordinate  variable  is  (k  /  x  e)  ^ 1  ^ 1  ~  *  The  abscissa  represents  t/tr. 
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We  now  consider  problems  for  steady  shocks  in  materials  with  strain  curves 


of  the  form  (5.5.2).  Our  results  are  given  in  terms  of  u,  the  scaled  strain  of  equa¬ 
tion  (4.6.5)  normalized  by  its  value  at  the  discontinuity: 

u(t)  =  u(t)/u0  -  s(t)/£0.  (5. 5. IS) 

After  we  derive  an  equaton  for  u,  we  present  graphs  of  numerical  approximations 
to  its  solutions  based  upon  the  scheme  presented  above  for  convolutions  with  an 
exponential  normalized  compliance.  We  then  modify  the  scheme  for  power-law 
normalized  compliances  and  present  the  graphs  of  approximations  for  shocks  in 
such  materials. 

To  obtain  the  governing  equation  for  t T,  we  cancel  f(e0)  =  \e0  from  both  sides 
of  equation  (5.3.2)  and  use  (5.5.2)  explicitly  for  /: 

k(e{t)y  =  (A  -  l)s(f)  +  A(X  '*  (5.5.19) 

c\ 

(Recall  A  =  ( U/U0 )  .)  In  terms  of  the  scaled  r  train,  this  is: 

(«(<))  7  =  (A  -  l)u(<)  +  A(x'*  «)(<)•  (5.5.20) 

We  note  that  this  equation  reduces  to  (5.5.4)  governing  acceleration  waves  at  A  =  1. 
For  A  >  1,  the  jump  in  the  shock  wave  has  the  value: 

u0  =  u  (0)  =  (A  -  1)~,  (5.5.21) 

since  =  0  implies  that  the  convolution  in  (5.5.20)  vanishes  at  t  =  0.  W  hen  we 

divide  (5.5.20)  by  v0,  we  obtain: 

(«(  t))1  =  «(<)+  (A7'  *  «)(0,  (5.5.22) 

w  here 

F(0=-^—  X(t)  (5.5.23) 

A  —  1 
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T  ho  c’orrospond  in  g  iteration  srh<»mi»  for  I\ 


'Tills  111  a  tl  <•  li  r  !>. 


«  *  +  l  (  t )  -  (  M  k  (  0  +  (  X  '  *  »  k  )  (  t) )  ' 1  '  ■  k  -  0.1.2.  A  1  ( 5  I  I  I 

We  again  choose  the  quasi-elastic  solution  to  start  our  upper  hound  sequences 

T0(  t)  =  uQ{  t)  -  (//  (  /)  +  \(  t))  1-1  ( r, ) 

It  is  easy  to  show  that  a  starter  for  a  lower  hound  sequence  is  the  Heaviside  step 
fn  n  c  t  ion  : 

«\>( '  )=«’;,(  0  =//  (  0  (5  5  2u  i 

I'' or  the  exponential  compliance-  in  equation  (5.5.1).  our  normalization  -  f  the 
solution  by  its  j  u  m  |>  value  requires  specification  of  the  parameter 

'-rrp'- 

For  any  choice  of  c  >  \e  the  normalized  wavespeed  being  considered  is: 

c 

^  =  - .  ( 5 .  5 .28) 

c  -  Vc 

In  what  follows,  we  have  chosen  c  =  1  to  illustrate  shock  behaviour  depending  on 
the  nonlinearity  q.  This  necessarily  restricts  the  compliance  to  have  an  equilibrium 
value  \  ,  <  1:  however,  this  is  not  a  severe  restriction.  Materials  for  which  \  .  <g;  1 
are  nearly  elastic;  with  small  values  of  the  graphs  in  Figures  5.5  5.8  repre¬ 
sent  numerical  approximations  to  shocks  at  wavespeeds  just  above  the  accelera¬ 
tion  wavespeed  in  nearly  elastic  materials  Furthermore,  values  of  approaching 
unity  represent  more  strongly  viscoelastic  solids  having  equilibrium  compliances 
near  J ,  =  2 J0,  according  to  (5  5.3).  For  such  materials,  the  graph-  represent  high 

8  2 


speed  cases,  >  3>  1 


plotted  in  each  curve,  rnr  slmcks.  the  <|  11  asi-eiast  ir  dilution  is  asymptotic  to  tin- 
solution  for  t  1  0  (§§5.1)  It  is  evident  from  these  figures  that  the  <1  II  asi-elasl  !■• 
solution  is  acceptable  as  an  approximation  for  a  quadratic  nonlinearity  with  this 
compliance;  it  is  very  good  for  q  5.  However,  for  q  =  f> /  1 .  it  is  no  better  than  a 
loose  upper  bound,  except  in  the  small-time  asymptotic  region 


Figure  5.8.  The  quasi-elastic  solutiou  and  the  superimposed  6-th  iterates  in  the 
upper  and  lower  bound  sequences  representing  the  approximate  numerical  solution 
for  the  strain  shock  wave  with  7  =  5,  an  exponential  compliance,  and  c  ~  1  T  !>  • 
ordinate  variable  is  e/e0.  The  abscissa  represents  t/tr. 


For  the  exponential  compliance  considered  thus  far,  we  could  have  success¬ 
fully  used  a  piecewise  linear  approximation  to  the  entire  integrand  in  the  convo¬ 
lution;  i . e . ,  the  familiar  trapezoid  quadrature  rule  [5.3],  We  chose  to  integrate  the 
iterates  exactly  with  d\  in  anticipation  of  the  singular  problems  which  we  now 
consider,  involving  the  power-law  normalized  compliance: 

X(t)=Xotp>  P  e  [0,1).  (5.5.29) 

The  numerical  scheme  is  modified  by  the  evaluation  of  the  integrals  / k  in  equa¬ 
tion  (5.5.16).  The  non-dimensional  time  in  this  case  is  xl^t.  We  make  the  change 
of  variables  x^J V t  — ►  t  in  (5.5.16)  and  use  d\(r)  =  prp_1  dr.  The  result  for  Ik  after 
simplification  is: 

h  =  [(l  ~  +  k  <pf_  t+1]  •  (xk  ~  X'jt-i) 

-  (ipf_*+1  -  V?-k)~ 7/,P  (fcP+1  “  (*  “  0P+1)-  (5.5.30) 

For  the  numerical  scheme  using  (5.5.30),  we  again  generate  upper  bound 

and  lower  bound  sequences  which  are  started,  respectively,  by  uq  and  L ^  tTjr_ . 
The  resulting  approximate  solutions  for  shocks  in  a  nonlinear  power-law  material 
are  illustrated  in  Figures  5.9-5.17.  These  figures  include  all  combinations  of  non¬ 
linearity  7  =  5/4,  2,  and  5  with  compliance  powers  p  =  1/10,  1/2,  and  9/10.  For 

these  figures,  we  used  401  points  in  the  computations,  all  of  which  are  plotted  for 

each  curve.  We  have  removed  the  wavespeed  from  the  the  problem  by  plotting  u 
against  the  non-dimensional  time  [ A x o / ( —  l)|1/,pf-  It  is  apparent  from  these  figures 
that,  for  the  problems  considered,  the  quasi-elastic  solution  is  better  as  an  approx¬ 
imation  to  the  solution  in  a  power-law  material  when  the  nonlinearity  7  is  large,  or 
when  the  power  p  is  small,  or  both. 
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Figure  5.9.  The  qu  asi-eiastic  solution  (v.p  curve)  and  the  su  pcim  posed  50-th 
iterates  in  the  upper  and  lower  bound  sequences  representing  the  approximate  nu¬ 
merical  solution  for  the  strain  shock  wave  with  7  =  5/4  and  p  =  1/10.  The  ordinate 


variable  is  e/c0.  The  abscissa  represents  [Ax„/(A  -  l)]'/p( 


Figure  5.10.  The  quasi-elastic  solution  (top  curve)  and  the  superimposed 
46-th  iterates  in  the  upper  and  lower  bound  sequences  representing  the  approxi¬ 
mate  numerical  solution  for  the  strain  shock  wave  with  ~i  =  5/4  and  p  =  1/2.  The 
ordinate  variable  is  e/e0.  The  abscissa  represents  [Ax0/(A  -  I)|1^p<- 


Figure  5.12.  The  quasi-elastic  solution  and  the  superimposed  H-th  iterates  in  the 
upper  and  lower  bound  sequences  representing  the  approximate  numerical  solution 
for  the  strain  shock  wave  with  7=2  and  p  =  1/10.  The  ordinate  variable  is  e/s0- 
The  abscissa  represents  [Ax<j/(A  -  l))l^pt. 


Figure  5.15.  The  quasi-elastic  solution  and  the  superimposed  6-th  iterates  in  the 


upper  and  lower  bound  sequences  representing  the  approximate  numerical  solution 
for  the  strain  shock  wave  with  7=5  and  p  =  1/10.  The  ordinate  variable  is  e/e0 
The  abscissa  represents  [Ax0/(A  —  l)]l/p<. 


95 


Figure  5.17.  The  quasi-elastic  solution  (top  curve)  and  the  superimposed  6-th 


a 


iterates  in  the  upper  and  lower  bound  sequences  representing  the  approximate  nu¬ 
merical  solution  for  the  strain  shock  wave  with  7=5  and  p  =  9/10.  The  ordinate 
variable  is  e/e0.  The  abscissa  represents  [Ax0/(A  -  l)]1/p<. 
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CHAPTER  6:  SOLUTIONS  BELOW  THE 
ACCELERATION  WAVESPEED 

6.1  Prelude. 

In  this  chapter,  we  consider  problems  for  non-negative  steady  strain 
waves,  e(t),  traveling  at  speeds  slower  than  that  of  the  acceleration  wave,  i.e.,  for 
U  <  UQ.  We  have  shown  in  Chapter  3  (§§3.3)  that  there  are  no  non-negative  non¬ 
trivial  steady  waves  which  travel  at  or  below  the  equilibrium  wavespeed  Ue.  From 
Chapter  4  (§§4.5),  we  know  that  when  U  £  (Ue,U0)  there  are  continuous  solu¬ 
tions  e(  <)  which  are  asymptotic  to  ert  as  t  — ►  —  oo  with  a  value  of  r  which  depends 
on  U  and  the  material  compliance.  In  section  6.2,  we  consider  the  problem  for  such 
a  solution  when  the  wavespeed  is  perturbed  from  the  limiting  value  Ue.  We  denote 
by  6  a  small  parameter  and  obtain  solutions  e(  t;  6)  which  are  asymptotically  correct 

rt  A 

as  6  [  0  when  U  =  t/e  ( 1  +  5).  In  section  6.3  we  produce  numerical  solutions  for  all 
wavespeeds  U  <  U0  in  power  law  materials  with  quadratic  nonlinearity. 

6.2  U2  =  U2(  1  +  5). 

We  first  derive  the  leading  term  in  an  asymptotic  series  for  e(t:6)  as  6  [  0 
when  the  strain  curve  /(e)  is  specified  no  more  completely  than  in  equation  (2. 3. IS). 
In  order  to  obtain  the  next  term  in  such  a  series,  or  to  just  estimate  the  size  of  the 
next  term,  we  need  to  further  specify  /.  We  do  so  in  the  rest  of  this  section  which 
we  devote  to  strain  curves  which  are  accurately  represented  by  a  power  series  for 
small  e.  For  these  materials,  we  obtain  the  first  two  terms  in  the  asymptotic  series 
for  e  along  with  an  estimate  of  the  size  of  the  third  term. 
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where  we  have  used  the  continuity  of  e  at  these  wavespeeds  (§§3.3)  in  the  last  line 
of  (6.2.7).  Thus,  the  governing  equation  (6.2.1)  becomes: 

f(i(0))  =  U2(J*i')(6).  (6.2.8) 

The  change  of  variables  6  —  rj  — ►  r;  in  the  integral  yields 

/oo 

e'(0  -  r,)j(r,)dV.  (6.2.9) 

o 

We  use  (5.1.8)  to  write 

J(9)  =  Je[ll(0)  -  j>(0)]  (6.2.10) 

in  equation  (6.2.9)  and  evaluate  the  first  integral  to  obtain: 

f  °° 

(J  *  i')(B)  -  Je£(0)  -  Je  e'(0  -  17 )  (*7 )  dtj.  (6.2.11) 

J  0 

We  assume  that  e  is  sufficiently  differentiable  to  do  the  following.  We  expand 
i' (6  —  t])  about  0  and  obtain: 

(J  *£')(6)  =  Je£(0)-Jej  fe'(^)-,e,/(<l)+%2e',V)-  +  •  1  tM  r,)dr,.  (6.2.12) 

0  L  21 

Let  us  define 

/oo 

Tj1>(r)dT-  j  =  0,1,2,  ...  ,  (6.2.13) 

0 

so  far  as  these  moments  exist.  Then, 

f  00  ■  . 

I  dV  =  r]+lrl>j.  (6.2.14) 

J  0 

Using  (6.2.13)  and  (6.2.14)  in  (6.2.12),  we  have: 

(J  *  e)(0)  =  Je  £(0)  -  ripo£'{0 )  +  r2ipii"(8)  -  —r3tp2e"'(0)  +  -  •  •  •  (6.2.15) 

If  the  moments  ip  •  fail  to  exist  beyond  some  value  of  j ,  the  terms  which  do  make 
sense  give  the  asymptotic  behaviour  near  r  =  0  (i.e.,  near  5  =  0)  (2.5).  In  what  fol¬ 
lows,  we  assume  that  the  moments  exist  so  far  as  we  need  them  for  computation. 


Wherever  the  small  parameter  r  appears  explicitly  in  (6.2.15),  we  need  to 
replace  it  by  an  expression  in  the  small  parameter  6.  Additionally,  we  will  ex¬ 
pand  £  in  a  series  in  6,  the  form  of  which  series  depends  on  the  strain  curve  /.  To 
first  determine  how  r  depends  upon  6,  we  use  equation  (6.2.4)  and  the  equalities 
U2  =  U2{1  +  6)  =  J~l(\  +  6)  to  write: 

r7(r)  =  Je(l  -  5  +  52  -  <53  +  •••).  (6.2.16) 

Recall  that  J  (r)  is  the  Laplace  transform  of  J(f).  We  transform  equation  (5.1.8) 
to  obtain : 

rJ(r)  =  Je[  1  -  r^(r)j.  (6.2.17) 


For  small  r: 


1>(r) 


-S 


e  rt  ip(  t)  dt 


n  L 


1  -  rt  +  —  t  -  +  •  ■  •  j  ip(t)  dt 


=  V'O  —  •01*'  +  — V»  2r  -  + 


(6.2.18) 


with  the  use  of  (6.2.13). 

Using  (6.2.18)  in  (6.2.17),  we  have: 


rJ(r)  =  Je 


1 


1  -  V'0r  +  lr“  -  ~V’2r  + 


(6.2.19) 


Comparison  of  (6.2.19)  with  (6.2.16)  suggests  a  power  series  expansion  of  r  in  6: 

r  =  6  r  j  4-  S2r2  +  63r3  +  -  •  •  (6.2.20) 

When  we  use  (6.2.20)  in  (6.2.19)  and  equate  the  result  with  (6.2.16),  we  find  that 


the  coefficients  of  the  powers  of  6  in  the  expansion  of  r  are: 

1 

rl  =  — 

V'O 


r2  =  rl(^lri  ~  l) 

r3  =  »T  (  -  +  2^lrlr2  +  • 


(6.2.21) 


and  so  on . 


3 


and  fj  is  determined  to  within  an  integration  constant  from: 


-  il  +  kij  =  0. 


(6.2  29) 


To  integrate  (6.2.29),  we  note  that  it  is  separable, 

1  d£  j 


1  -  kej  1  e‘l 


dd 


(6.2.30) 


and  may  be  put  in  the  form 


dr) 


1  -  he”/0 


=  de. 


(6.2.31 ) 


where  a  is  given  in  (6.2.27)  and 

i)  -  In  £  i . 

Upon  integration  of  (6.2.31),  we  have: 

ke’1'a  = 


(6.2.32) 


t(»-»  o)/« 


(6.2.33) 


]  4.  e[8- 0o)/°  ’ 

w'here  6 0  is  the  integration  constant.  In  this  problem,  the  time  origin  is  arbitrary. 
We  choose  to  set  this  origin  with  0O  =  0.  Using  (6.2.27)  and  (6.2.32)  to  write  the 
result  in  terms  of  fj,  we  have: 

1  eb-1)*  } 


*l(*) 


fl  eb-W  1 

"  U  ’  1  +  I 


+ 

This  may  be  written  in  terms  of  the  hyperbolic  tangent  as: 


(6.2.34) 


it («)  = 


2k 


1  +  tanh 


0  -  1 


_i _ 

i-i 


(6.2.35) 


With  the  use  of  (6.2.5),  we  observe  from  (6.2.34)  that  ej  has  the  asymptotic  form 
required  by  (6.2.3),  to  within  the  arbitrary  time  shift  implicit  in  the  multiplicative 
constant  k~  1 


We  now  consider  this  problem  in  more  detail  for  strain  curves  amenable  to  a 
power  series  expansion: 


■  1  Z  Ml . *  . .  S  I  ■  1  1  I  I  I  11  i  i  I 


4  a  4*a  <  a  i*a  >  a. | '»  4»«.t  >>  *■*  j.|jj  | 


i.»  a„»  i.<  t-*  iw»  *-■ 


In  this  case,  it  is  reasonable  to  assume,  at  the  outset,  that 


i(»)  =  Si  i(6)  +  62£2(8)  +  0(i3), 


(6.2.37) 


where  the  functions  £\(9)  and  €2(6)  and  their  derivatives  are  0(1)  with  respect  to  8 . 


We  use  (6.2.37)  in  (6.2.36)  and  obtain: 


/(e)  -  i  =  <52ae2  +  53(2ae,e2  +  /5e3)  +  0(<54). 


(6.2.38) 


Similarly,  from  (6.2.37)  and  (6.2.22)  we  have: 


U2(J  *  £')  -  £  =  82(£l  -  i\) 


+  (53(e2  -  e2  +  riV'i(^j  -  «i)) 


+  0(«4). 


(6.2.39) 


We  equate  (6.2.38)  and  (6.2.39)  according  to  (6.2.8)  and  obtain  equations  which  ej 


and  i2  must  satisfy.  From  the  0(6 2)  terms,  we  get  equation  (6.2.29)  for  this  case: 


-/  .  .  .2  „ 

el  ~  el  +  a£i  — 


(6.2  40) 


The  0(63)  terms  produce  an  inhomogeneous  differential  equation  for  e2  with  vari¬ 


able  coefficients  which  are  continuous  functions  of  8: 


£'2(9)  -  a(6)£2(9)  =  6(0), 


(6.2.41) 


w  here 


a(9)  =  (1  -  2al1(0)) 


b(8)  =  TjV'i (£i(9)  -  4(0))  -  />e-*(0). 


(6.2.42) 


We  obtain  £t  immediately  from  equation  (6.2.35): 


«l(0)  =  ~[l  +  tanh(0/2)j . 


(6.2.43) 


The  solution  for  e2  is  the  sum  of  a  particular  integral  and  the  homogeneous  solution: 


il(0)  =  4(9)  +  4(°)- 


(6.2.44) 


The  solution  for  the  homogeneous  equation  is  in  general: 


*£(*)  =  C  exp  /  a(Z)dZ  >, 


(6.2.45) 


where  C  is  the  integration  constant.  With  a(6)  and  £i(0)  given  by  (6.2.42)  and 
(6.2.43),  the  integration  is  for  >s  straightforward.  The  result  is: 

e^O)  =  C  cosh' 2 (6/ 2).  (6.2.46) 


For  eAQ),  we  have  in  general: 


W  =  /  exP  |  /  «(t)  dZ  |  M*?)  drl- 


(6.2.47; 


Using  £i(0)  in  equations  (6.2.42)  for  a{0)  and  6(0),  we  have  after  much  simplifica¬ 
tion,  the  intermediate  result: 

3  r  6/2 

£V2(6)  —  cosh  2(0/ 2)  X/  cn  J  sinhn(£ )  cosh2 ~ ” ( £)d£,  (6.2.48) 


w  here 


*1  , 

0  \ 

1  / 

'^1 

3  p 

~ ,  + 
^0 

2a3  /  ’ 

c2  =  - l 

2a  ' 

2a3 

*1 

l  1 

f  V»i 

0 

*1 

2o3  / 

C3  =  57 1 

v  ^0 " 

2a3 

(6.2.49) 


for  which  we  have  used  (6.2.21)  to  write  the  c„  in  terms  of  the  compliance  mo¬ 
ments,  rj) 0  and  V'l.  defined  in  (6.2.13).  Again,  the  integrations  in  (6.2.48)  are 
straightforward.  After  integrating  and  simplifying  the  result  for  e^,  we  add  it  to  £ 
according  to  (6.2.44).  The  result  for  ^2!®)  >s: 

e2(0)  =  —  — -j  tanh  {0/2)  4-  tanh‘(0/2)j 

+  {^(^t}{C+  ii(^^)((</2)+  [«=o»h(*/a)])}.,«.2.so| 


To  verify  that  £2(6)  has  the  correct  asymptotic  property  for  0 


—  00  ,  we  use 


the  following  identities: 


cosh  x  = 


ex  +  e  1 


(6.2.51) 


tanh  x  ~ - 1 . 

ex  +  e~  z 


For  large  negative  x: 


cosh  x 


x  — ►  —  00  . 


(6.2.52) 


tanh  x  ~  2e  1  -  1, 


Thus, 


tanh (6>/2)  +  tan h 2 (0 /2 )  ~  -2e\ 


cosh~2(0/2)  ~  4efi 


8  — ►  —  00  . 


(6.2.53) 


In  [  cosh (0/2)]  ~  - (6 / 2)  -  In  2, 


Using  (6.2.52)  and  (6.2.53)  in  (6.2.50),  we  obtain: 


00.  (6.2.54) 


As  in  the  solution  for  £1(8),  the  integration  constant  C  serves  to  define  the  arbitrary 


time  shift  of  the  solution. 


With  and  £2  given  in  (6.2.43)  and  (6.2.50),  our  solution  for  £{6)  is  com¬ 


plete,  so  far  as  we  have  expanded  it  in  equation  (6.2.37).  Since  £j  and  £2  are  bounded 


in  8,  the  expansion  is  uniformly  asymptotic.  To  go  on  to  higher  orders  in  the  so¬ 


lution,  one  merely  continues  the  expansion  of  (6.2.37)  with  terms  6C  £ 3(0)  +  0(6*) 


and  obtains  an  equation  for  £3(0)  in  terms  of  £j  and  £2- 


As  a  final  remark  on  the  solution  we  have  obtained,  we  show  that  it  produces 


the  correct  equilibrium  value,  to  within  0(63).  From  (6.2.43)  and  (6.2.50),  respec¬ 


tively,  we  have: 


lim  £  j  (0 )  =  -I — 
$  —  +  00  ct 


(6.2.55) 


lim  c2(0)  = - r. 


S-*  +  00 


Thus,  the  equilibrium  value  for  e(0)  from  (6.2.37)  is: 


£1 


6  , 
it  =  -  -  Hr  +  0(6 3). 


(6.2.56) 


We  recall  from  section  3.3,  equation  (3.3.2),  that  the  equilibrium  value  satisfies: 


/(*.)  =  V*Jeit 


(6.2.57) 


Using  U2Je  =1  +  6  and  the  expansion  of  /  in  (6.2.36),  we  have  from  (6.2.57): 


aie  +  Pi]  +  0(i])  =  6. 


(6.2.58) 


When  we  solve  (6.2.58)  for  ee  in  terms  of  6,  we  obtain  equation  (6.2.56),  the  required 
result. 


6.3  Power  Law  Materials  with  Quadratic  Nonlinearity. 


We  now  obtain  steady  wave  solutions  at  all  wavespeeds  below  that  of  the 
acceleration  wave  for  materials  having  a  purely  quadratic  nonlinearity  in  the  strain 


/ (e)  =  e  +  are  ,  a  >  0. 


(6-3.1) 


Our  solutions  apply  for  two  different  classes  of  material  compliance  models.  The 
first  class  is  characterized  by  compliances  J  which  have  a  non-zero  initial  jump,  J0, 
followed  by  a  power  law  norm  alizcd  compliance  (§§4.6): 


J(t)  =  J0[H(t)  +  CP(t)\ 


pm  = 

p! 


P  €  (0,1], 


(6.3.2) 


where  p!  is  defined  in  equation  (2.4.11)  and  C  is  a  positive  constant.  Another  class 


of  material  models  is  defined  by  a  power  law  compliance  with  no  initial  jump: 


J(t)  =  cp(t),  p  e  (o,  i], 


(6.3.3) 


V.i  VtVtV  r. 


i  ivr,  iV»V*' 


S*  V?.« 


1  *i» 


where  C  is  a  positive  constant.  In  such  a  model,  p  is  the  log-log  slope  of  the 
time  dependence  of  J  and  is  fairly  constant  over  broad  ranges  of  t  in  real  materials 
(see,  e.g.,  Ferry  [2.3]).  In  either  model,  p  =  0  corresponds  to  purely  elastic  ma¬ 
terials.  Additionally,  neither  compliance  model  has  a  finite  equilibrium  value;  the 

O 

equilibrium  wavespeed  =  1  /Je  is  therefore  zero.  Steady  waves  can  exist  in  these 
materials  for  all  non-zero  values  of  wavespeed.  Since  the  power  law  material  in 
equation  (6.3.3)  has  no  initial  jump,  it  cannot  support  steady  shock  or  acceleration 
wave  solutions;  the  acceleration  wavespeed  U 0  =  1  /J0  is  infinite  for  J0  —  0.  For 
such  materials,  the  steady  wave  solutions  produced  in  this  section  apply  for  all  finite 
wavespeeds  U  €  (0,oo). 

For  either  compliance,  we  reduce  the  problem  to  that  of  solving  the  same 
nonlinear  integral  equation  with  a  singular  kernel.  For  the  normalized  power  law 
model  (6.3.2)  with  the  strain  curve  of  equation  (6.3.1),  the  governing  equa¬ 
tion  (3.1.10)  yields: 

(1  -  A)e  +  as2  =  \C(P  '  *  e)(  t),  (6.3.4) 

for  which  we  recall  the  definition  of  A  given  in  equation  (5.3.3): 

A  =  U2J0  =  (U/U0)2.  (6.3.5) 

For  the  wavespeeds  of  this  section,  we  have  A  <  1.  We  introduce  the  scaled  strain 

u(t)  =  (6.3.6) 

1  A 

and  obtain  from  (6.3.4): 

A  C 

«(<)  +  «2(<)  =  - -(/>'*  «)(<)  (6.3.7) 

1  —  A 

The  positive  constant  multiplying  the  convolution  may  be  absorbed  by  introducing 


the  scaled  non-dimensional  time 


We  then  have 

«’('?)  +  »2(n)  =  {P  '  *  «)(v),  (6.3.9) 

w  here 

u(t])  =  u(  t).  (6.3.10) 

Similarly,  for  a  pure  power  law  compliance  defined  in  (6.3.3),  equation  (6.3.9)  ob¬ 
tains  when  u  and  t)  are  defined  as 

u(t)=ae(t)  and  »/  =  (U2C) 1/p  t.  (6.3.11) 


For  both  classes  of  materials,  equation  (6.3.9)  is  the  nonlinear  integral  equation  to 
be  solved  for  u.  The  kernel  P'(i})  is  singular  for  p  €  (0,  1).  If  p  =  1,  the  kernel  is 
the  Heaviside  step  function,  with  which  convolution  is  a  pure  integration.  In  this 
case,  we  have 


[  n 

«(>?)  +  “2(»?)  =  J  «(»)  ds • 

~  OO 

For  any  differentiable  nonlinear  function  /(«),  the  integral  equation 


(6.3.12) 


/( 


.»./*  • 


(T)  dT 


(6.3.13) 


may  be  differentiated  to  obtain 


/,(u)d(ln  u)  =  dr. 


(6.3.14) 


Integration  of  (6.3.14)  gives  u  implicitly: 


J  f'(s)  {ds/s)  =  t  -  t0,  (6.3.15) 


for  an  arbitrary  shift  t0.  For  our  quadratically  nonlinear  problem  with  p  =  1,  the 


exact  solution  from  (6.3.12)  is  therefore: 


w  ’vi 


w  here 


K  i 


(*  +  l)p’ 


/r  =  0,1,2,  ...  . 


(6.3.23) 


With  the  observation  that  K  q  =  1,  we  see  that  c0  is  arbitrary;  we  choose  cq  =  1  to 
obtain  the  solution  which,  according  to  (6.3.18),  is  asymptotic  to  en  for  rj  — *  —  oo  . 
Since  {e*l,}£°=i  's  a  linearly  independent  set,  the  coefficients  in  the  series  of  (6.3.18) 
are  determined  from  the  recursion  relation: 


CQ  —  1 


ck  = 


~(c  *  c)k_  i 

1  -  K , 


A-  =  1,2,3,  ...  . 


(6.3.24) 


We  note  that  (c*c)o  =  co=  1-  Furthermore,  A' k  6  (0,  1)  for  k  >  1.  Thus,  Cj  <  0.  It 
follows  by  induction  that  the  coefficients  ck  alternate  in  sign,  since  the  summation 
(the  discrete  convolution)  involves  terms  which  are  either  all  negative  or  all  positive. 


Given  the  series  for  tT  defined  by  (6.3.18)  and  (6.3.24),  we  must  show  that 
it  converges  absolutely  on  a  non-empty  interval  (-oo,»/j-].  Viewing  (6.3.18)  as  a 
power  series  in  e’1 ,  we  are,  in  effect,  looking  for  a  non-zero  lower  bound  on  the  radius 
of  convergence  of  the  series 

oo 

Si(t)  =  52  \ck\zk.  (6.3.25) 

k=  0 

It  is  difficult  to  directly  estimate  the  growth  of  the  coefficients  in  this  series  We 
note  for  later  comparison  that  the  coefficients  obey 


lcfcl'>  (c*c)/t— i>  A  —  1,2,... 


We  instead  consider  the  series 


00 

S2(  z)  =  E  bkzkt 
k=  0 

with  strictly  positive  coefficients  defined  by 


bk  =  (-l)*c*(l  -  Ki)k,  k  =  0,1,2,  ... 


(63  26) 


(6.3  27) 


(6.3.28) 


If  52  converges  for  |z|  <  R,  then  Si  does  so  for  |z|  <  (1  —  Ki)R.  Consequently,  the 
series  in  (6.3.18)  for  u  converges  absolutely  for  t]  <  In  [(1  -  Kk)R],  To  show  the 
convergence  of  S 2,  we  first  determine  the  recursion  relation  for  its  coefficients.  We 
will  then  show  that  these  coefficients  are  dominated  by  those  of  a  power  series  with 
a  non-zero  radius  of  convergence.  Using  equations  (6.3.24),  (6.3.19),  and  (6.3.28), 

we  obtain  the  recursion  relation: 

6q  =  1 

1  -  Kj 

bk  = - —  (6  *  6)*_i.  A*  =1,2,....  (6.3.29) 

1  ~  Kk 

Observe  that  K  k  <  K\  for  all  k  >  2.  Therefore, 

bk  <  (6  *  b)k_  1,  k  =  1,2,...,  (6.3.30) 

with  equality  only  at  k  =  1.  We  have  transformed  series  whose  coefficients 
are  bounded  below  by  their  convolution,  equation  (6.3.26),  into  series  S 2  for  which 
the  coefficient  convolution  provides  an  upper  bound.  Consider  now  the  se¬ 
quence  {dfc}“_0  of  positive  numbers  defined  by: 

d0  =  1 

(6.3.31) 

dk  =  (d  *  d)k_l. 

Since  60  =  d0  =  1,  it  follows  from  (6.3.30)  and  (6.3.31)  that 

bk  <  dk,  k  =  0,1,2,  ...  ,  (6.3.32) 

with  equality  only  at  k  =  0,1.  The  sequence  {dk}<£’_Q  dominates  the  se¬ 

quence  {6jt}£°_0-  We  now  exhibit  a  function  which  generates  a  power  series  having 
a  non-zero  radius  of  convergence  with  coefficients  which  obey  the  recursion  rela- 


with  the  branch  of  the  square  root  which  is  positive  when  its  argument  is  real  and 

positive.  For  \z\  <  1/4,  we  may  write  g  as  a  convergent  power  series: 

oo 

g[  z)  =  E  Ok2* .  (6.3.34) 

k  =  0 

where  the  sequence  of  coefficients  {<7fc}£°_0  can  be  determined,  for  example,  by 

expanding  the  square  root  in  a  binomial  series.  In  particular, 

90  =  1 


and 


9/f  >  0  for  all  k. 


(6.3.35) 


The  series  is  absolutely  convergent  and  we  may  square  it: 

00 

g2( z)  =  E  (g  *  g)kz>c-  (6.3.36) 

k= o 

From  (6.3.33),  we  eliminate  the  radical  to  obtain 

z(g  —  l)  =  z2  g2  ■  (6.3.37) 

Using  equations  (6.3.34)  and  (6.3.36)  in  (6.3.37),  we  have 

00  oo 

E  9kzk+l  =  E  (9*  g)k*k  +  2>  (6.3.38) 

*=1  *  =  0 

where  we  have  used  jo2°  =  1-  Upon  adjustment  of  the  indicies  on  the  right  sum- 
m  ation  ,  we  see  that 

9k  =  (9  *  g)k-  l>  k  =  1,2,  ...  ,  (6.3.39) 

in  view  of  the  linear  independence  of  {zfc}£°_ ,.  Since  g  q  =  1,  the  coefficient 

sequence  ‘s  identical  to  the  sequence  {djt}£°_0.  Thus,  the  convergence 

for  |  z  |  <  1/4  of  the  series  for  g  which  dominates  implies  the  convergence  of  So 
in  the  same  circle.  Hence,  for  any  which  satisfies 


<  In  [(1  —  Ki)/A] , 


(6.3.40) 


the  series  for  tT( r; )  in  (6.3.18)  converges  absolutely,  at  least  for  all  t}  <  t)j 
Since  1  —  K  \  <  1,  we  have  t)p  <  0.  We  note  that  the  series  S  j  diverges  for  |;|  >  1/4. 
For,  equations  (6.3.26)  and  (6.3.39)  with  go  =  1  imply  that  the  coefficients  of  Si(  z) 
dominate  those  of  the  series  for  g(  z),  which  series  for  g  diverges  outside  the  cir¬ 
cle  |r|  =  1/4.  While  the  series  for  tT  may  converge  absolutely  for  some  r)  larger 
than  tjf,  it  diverges  absolutely  for  r)  >  In  ( 1  /  4 ) ,  for  all  values  of  p. 


We  have  thus  far  solved  the  problem  on  the  half-line  t]  <  1/7-,  where  «(»/)  is 
given  exactly  by  the  alternating  series  defined  by  equations  (6.3.18)  and  (6.3.24).  We 
construct  the  solution  for  i)  >  rjT  by  iteration  of  monotone  sequences  of  monotone 
functions,  much  as  we  did  for  the  earlier  constructions  of  shocks  and  acceleration 
waves.  For  the  current  problem,  however,  each  iterate  involves,  for  »/  >  i]p.  a 
contribution  from  the  series  solution  for  rj  <  rjp.  The  iteration  scheme  is: 

«*+l(»j)  =  F((p  '  *  «*)('/))  -  (6.3.4  1 ) 


_  _ O 

where  we  now  denote  by  F  the  inverse  of  «  +  w  .  We  first  show  that  an  increasing 


(lower  bound)  sequence  is  started  by 

«o('/) 


=  /  v_(ri),  tl  <  t]T\ 

l  «(»/r).  1  ^  »?r; 


(6.3.42) 

where  u  is  the  known  series  solution  for  i)  <  rjp.  Recall  that  u  is  monotone,  accord¬ 
ing  to  Theorem  4.5.1.  For  t;  >  r/p, 


( P  '  *  « 0 )  ( '/ )  =  /  jP  (*/  -  a)  du0(s) 
*  -  00 
f  nr 

=  J  P(r)  -  s)  du0(s) 


■  OO 

/  — \ 


>(/>'*«)(  t,T),  (6.3.43) 

where  we  have  used  the  increasing  nature  of  u0  and  P  along  with  the  fact 


that  du0  =  0  for  >/  >  i)p.  Since  F,  too,  is  increasing,  we  have 


ui  (tj)  =  F  (  (P  '  *tic)(i|)j  >  F  (  (P  '  *u0)(Vr))  =  =  uo(*l),  */  >  V  T  (6.3.44) 


and  the  assertion  is  proven.  Furthermore,  u  [  is  bounded  above  by  e,  in  view  of 
Lemma  4.5.1.  Theorem  4.5.1  ensu res  that  iteration  of  tT0  converges  to  the  solution  t7. 

Similarly,  a  decreasing  (upper  bound)  sequence  is  started  by  the  discontinuous 
starter  defined  by: 

*•(’)  =  {  \(,)'  ”  ;  (6.3  15) 

1  e  ,  i]  >  ii  t  ■ 

For  r)  >  t)T , 

[(P  '  *  )  ( »; )  ^  <  P'  *  e’}  —  ^  (6.3.46) 

according  to  Lemma  4.5.1,  since  ii  (t /)  <  en .  Thus, 

«l(»?)  <  P{en)  <  en  =  tT0('/),  V  >  nr-  (6.3.47) 


All  iterates  in  the  resulting  upper  bound  sequence  are  bounded  below  by  the 
lower  bound  starter  of  equation  (6.3.42),  since  ITj  ( ty )  is  so  bounded,  in  view  of 
Lemma  4.3.4.  Accordingly,  the  upper  bound  sequence  converges  to  the  solution  u. 
We  remark  that  the  discontinuity  in  the  upper  bound  starter  is  not  present  in  the 
succeeding  iterates,  since  P  is  continuous. 

We  again  resort  to  a  simple  numerical  scheme  to  generate  graphs  of  approxi¬ 
mations  to  the  upper  and  lower  bound  sequences  discussed  above.  The  results  we 
present  are  not  for  solutions  asymptotic  to  e’1  for  t/  —*  —  oo  .  Rather,  we  consider 
th  e  sh  ifted  solu  tions 

«(»/)  =  “  I'/rl).  (6.3.48) 


so  that  the  series  solution  for  u  is  always  used  for  7]  <  0.  We  have 

00 

«(»/)=  E  cke{k+l)\  »/  <  0,  (6.3.49) 

t  =  0 

w  here 


-  (*+ 1  )|»?tI 


(6.3.50) 


and  t]f  is  chosen  according  to  (6.3.40).  Any  iterate  in  the  sequence  of  bounding 


functions  has  the  form: 


(6.3.51) 


<p{ i,)  =  J  «('/),  V  <  0; 

(  >p(n),  n  >  0. 

The  convolution  with  P1  is: 

[  00 

J  <p(>i  -  s)P'{s)  ds 
0 

/r]  r  oo 

<p(i)  —  s)P' (s)  d s  +  J  u(t]  —  s)P'(s)  ds.  (6.3.52) 


From  equations  (6.3.2)  and  (6.3.49),  we  see  that  the  last  integral  in  (6.3.52)  is  a 
sum  of  integrals  of  the  form: 


I(ri) 


'-j 


V-  n 

With  the  change  of  variables  rs  —  £,  we  have 

-perq  r  oo 


er[,,-l]pSp~'  ds. 


(6.3.53) 


r~pern  f 

I(1)=  : - —  j  e-'e-'di. 

(p  -  1)!  ru 


(6.3.54) 


The  integral  in  (6.3.54)  is  an  incomplete  gamma  Junction.  In  the  notation  of 
Abramowitz  and  Stegun  (5.1|,  we  have 

r (p,  rp) 


/(»?)  =  r"per" 


r  (p) 


(6.3.55) 


for  which  the  connection  with  our  notation  in  (2.4.11)  is  provided  by 


p!  =  r  (p  +  1). 


(6.3.56) 


Using  equations  (6.3.49),  (6.3.53),  and  (6.3.55)  in  (6.3.52),  the  complete  expression 
for  the  convolution  for  p  >  0  is: 


,  1  ^  e(t+1)rj 

(P  *  p)(v)  2^  ck  -~--pr(p,(fc  +  l)v) 

'  *  =  0 


'-r 


p(i;  -  s)  dP(s). 


(6.3.57; 
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The  alternating  series  in  this  result  can  he  approximated  arbitrarily  closely  by  a 
finite  number  of  terms  with  the  magnitude  of  the  error  bounded  by  that  of  the  first 
neglected  term.  To  approximate  the  incomplete  gamma  function,  we  use  the  con¬ 
vergent  series,  continued  fraction,  and  asymptotic  series  representations  given,  re¬ 
spectively,  in  equations  (6.5.29),  (6.5.31),  and  (6.5.32)  of  [5.1].  We  approximate 
the  convolution  integral  with  the  scheme  presented  in  equations  (5.5.15).  (5.5.16), 
and  (5.5.30)  of  this  work. 

Figure  6.1  illustrates  truncated  upper  and  lower  bound  sequences  for  a  quad- 
ratically  nonlinear  material  with  p  =  1/1.  Here  wc  chose  t/y  =  —  1  and  we  have  plot¬ 
ted  u(ij)  =  u"(»;  —  1)  against  i) .  It  shows  the  lower  bound  starter  of  equation  (6. 3. -12) 
and  the  first  19  iterates  of  the  sequence  it  generates,  along  with  the  upper  bound 
starter  of  equatio n  (6.3.45)  and  the  first  17  iterates  of  its  sequence.  For  each  approx- 
imate  iterate,  100  points  were  computed  and  plotted  for  »/  >  0.  For  all  figures  in 
this  section,  50  points  are  plotted  for  i/  <  0  where  the  alternating  series  solution  is 
shown.  Figures  6. 3-6. 6,  discussed  below,  were  prepared  with  200  points  computed 
and  plotted  for  ij  >  0.  As  in  the  graphs  of  Chapter  5,  the  large  space  between  the 
truncated  sequences  is  the  numerical  approximation  to  the  bound  within  which  the 
solution  is  to  be  found. 

We  expect  that  a  starting  function  which  is  similar  to  the  expected  form  of 
the  solution  will  yield  better  bounds  upon  iteration.  We  now  present  such  a  starter 
for  an  upper  bound  sequence.  We  define 

B  i  ( r/ )  =  c'1 ,  (6.3.58) 

and  Bjl'l)  as  the  solution  to: 

/?|(»/)  =  {r1  *  n2)(n)-  (6.3.59) 


The  solution  for  B 2  is 


B2(v)  = 


(2  p) 


(’I  ~  ’I 0 ) ? 


We  show  that  the  starting  function  defined  by 

«<>('/) 


=  /  #l(»/).  V  <  ’ll ; 
(  Bz(’l),  11  >  rj , , 


w  ith 


n  1  =  In  (—yy)  and  10  =  ’ll  - 


(6.3.60) 


(6.3.61 ) 


(6.3.62) 


starts  a  decreasing  sequence  upon  iteration.  These  values  of  r/()  and  r/  ]  produce  a 
continuously  differentiable  starter  u0.  From  the  iferation  scheme  (6.3.4  1).  we  have 


_  _ n  f  _ 

V  j  +  U  J  =  P  It  u, 


(6.3.63) 


_  _ o  _  _ n  _  _ 

If,  for  all  rj ,  iteration  produces  ti(  4-  u  “  <  u0  +  tt ~ ,  then  vt  <  1/  0  and 
a  decreasing  sequence.  It  is  therefore  sufficient  to  show  that  P  '  *  u0  < 
Since  P1  *en  —  eri ,  the  result  is  immediate  for  »/  <  »;  x .  For  i]  >  r/j, 


(P'*u0)(,,)  =  / 

« o ( s ) p '( n  -  s)  ds 

-  OO 

f  n' 

,  f " 

Bn(s)P'(tl  -  s)  ds 

=  J 

Bt  (s)P'(i,  -  s)ds  +  / 

-  OO 

‘  ’ti 

[  ’ll 

B  n  ( s )  P '  ( r/  -  ,< )  rf  a 

=  1 

Bl(s)P'(>l  -  s)  da  +  J 

-  00 

'to 

f  "■ 

-  J  Bn(s)P'(,l 

—  s  )  d  s . 

u'o 


starts 

_ O 

+  u 


0  ■ 


(6.3.64) 


'to 

Since  Bn  vanishes  for  r;  <  r/j,  we  can  extend  the  lower  limit  of  integration  in  the 

last  two  integrals  above  to  —  00  .  Thus, 

f  ’i  1  [1 

(P  '  *  u0)in)  =  J  («,  -  B2)(s)1>'(i,  -  *)ds  +  J  Bn{s)P'(,,  -  s)ds.  (6.3.65) 

-  OO  “00 

The  last  in  tegral  is  (P1  *  B^)  {’I )  =  B^(i]),  according  to  equ  at  ion  (6.3.59) .  T  he  re  fore, 


(P  '  *  u0)(l)  = 


Bn(l/)  + 


Bn  )(S)P'(’I  -  S)  ds 


< 


Bl(’l)  + 


f  n' 

J  {Bi  -  B2)(s)r'(,h  -  9)  da, 


(6.3.66) 


since  (Z?i  —  Zfo)(f/)  >  0  for  all  r;  /  i/ j  and  P'(’l)  is  a  decreasing  function.  The  last 


integral  in  (G.3.66)  is  P1  *  (B  j  -  Bn)[i/l)  =  B  j(Vl)  —  B~(i) ,).  Therefore, 

(p  '  *  «„)(»,)  <  Bl(,,)  +  «,(»,,)  -  «■;(.,,) 

<  /?$(r;)  +  W,('7i) 

=  P  a  ( t; )  +  P  2  ( ’7 1  ) .  ( G  .3  .6  7  ) 

since  u0  is  continuous  at  t;  j.  Now  /#  ._>  is  increasing,  Hence,  we  obtain: 

(P  '  *  «„)(./)  <  «.;;<»,)  +  fln(V) 

=  «”a(»/)  +  «o('/).  (G.3.GS) 


which  proves  the  assertion. 

With  this  result,  it  is  easy  to  show  that  an  upper  bound  starter  useful  for  our 
numerical  work  is  provided  by  the  discontinuous  function  defined  as: 

(  »"('/).  >7  <  '77-: 

«o(»7)  -  j  B,(t;)t  ijT  <  v  <  Vi  ( G  3 .GO ) 

1  yj(’7),  >7  >  *7 1 . 

with  »/  o  and  »/ 1  given  in  (6.3.62)  and  the  series  for  u  in  (6.3.18).  Figures  G  .2  G.G 
illustrate  approximate  bounds  based  upon  this  upper  bound  starter.  In  Figure  6.2 

(j>  =  1/4)  with  100  points  computed  and  plotted  for  r/  >  0,  we  see  that  the  early 

iterates  are  much  better  upper  bounds  than  those  started  by  the  discontinuous 
exponential  starter  defined  in  equation  (6.3.15)  and  illustrated  in  Figure  6.1 

Recall,  we  have  the  exact  solution  for  u  ( ;/ )  for  all  t/  in  the  single  case  when 

]>  —  1  It  is  given  in  equation  (6  3  16).  In  Figure  6.3.  w e  compare  the  results 

of  20  iterations  in  both  lower  and  upper  bound  sequences  with  the  exact  result 
to  provide  some  verification  of  our  numerical  scheme.  The  approximate  iterated 
solution  and  the  exact  solution  are  indistinguishable.  Figures  6.4  6.6  illustrate  our 


approximations  to  bounds  on  the  solutions  for  p  =■  1/2,  1/4,  and  1/10. 


it wvwwww  wjwi  wvu  i  ruwu  wj  hv  i*.  ■  w 


Figure  G.3.  The  exact  solution  for  p  =  1  and  the  superimposed  20-th  iterates  in 
the  upper  and  lower  bound  sequences.  The  ordinate  variable  is  ti  ( r/ )  =  u(t]  —  2). 
The  abscissa  represents  the  non-dimensional  time  p. 
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